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Tliis  thesis  presents  a  Time  Dei)enclent  Adaptive  Filter  (TDAF)  which  exploits 
the  cyclostationarity  of  digitally  modulated  communications  signals  and  seeks  to  im¬ 
prove  the  Signal  to  Interference  Ratio  (SIR)  and  Signal  to  Noise  Ratio  (SNR)  of  such 
signals.  The  TDAF  is  imbedded  in  a  computer  simulation  of  a  simple  communica¬ 
tion  system  consisting  of  a  data  source,  data  formatter,  pulse  shaping  filter,  BPSK 
modulator,  and  demodulator.  In  the  simulation  the  TD.AF  and  a  Time  Independent 
.Adaptive  Filter  (Tl.AF)  attempt  to  extract  the  Signal  of  Interest  (SOI)  from  noise 
or  interference.  The  criteria  of  Mean  Squared  Error  (MSE)  is  used  as  the  primary 
means  to  compare  the  performance  of  the  two  adaptive  filters.  Plots  of  MSE  im¬ 
provement  in  interference  and  MSE  improvement  in  noise  are  presented.  For  the 
case  of  interference,  the  improvement  is  measured  as  a  function  of  the  baud  rate  of 
the  intereference  signal,  and  carrier  frequency  of  the  interference  signal.  It  is  shown 
that  with  respect  to  the  TIAF,  the  TDAF  provides  up  to  12  dB  of  improvement. 
Bit  Error  Rates  (BER)  for  several  simulations  are  presented.  The  data  indicate  that 
signficant  improvements  in  BER  might  be  expected  when  a  TDAF  is  used  in  lieu  of 


a  TIAF. 
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FOR  COCHANNEL  INTERFERENCE 
REDUCTION 


/.  Introduction 


1.1  Background 

The  advent  of  frequency  reuse,  particularly  in  satellite  communications,  has 
made  the  task  of  recovering  digitally  modulated  signals  more  challenging.  Unoccu¬ 
pied  frequencies  in  the  electromagnetic  spectrum  either  do  not  exist  or  are  imprac¬ 
tical  for  use  due  to  their  high  frequency  or  poor  channel  performance  (in  the  O2 
and  H2O  absorption  bands,  for  instance).  A  single  communication  system  may  re¬ 
use  frequencies  up  to  six  times  through  polarization  and  spatial  (antenna  pointing) 
reuse  techniques.  INTELSAT  VI  is  an  example  of  such  a  system  (3:99). 

The  chissical  technique  of  applying  the  signal  of  interest  (SOI)  to  a  relatively 
narrow  passband  filter  can  be  of  limited  effectiveness  when  the  SOI  has  been  cor¬ 
rupted  by  another  signal  (or  signals)  not  of  interest  (SNOl)  whose  spectral  com¬ 
ponents  overlap  those  of  the  SOI.  Of  course,  with  only  knowledge  of  the  center 
frequency  and  spectral  width  of  the  SOI,  a  digital  filter  can  be  readily  designed  that 
limits  power  outside  the  spectral  region  of  interest  (.5:403-489).  The  impulse  response 
of  such  a  filter  is  fixed.  Hence,  if  the  characteristics  of  the  SNOI  changes,  the  filler 
cannot  appropriately  adapt  to  the  new  conditions. 

If,  on  the  other  hand,  the  filler  coefficients  arc  allowed  to  vary  so  as  to  minimize 
some  error  criteria,  the  impulse  response  will  no  longer  be  fi.xcd.  Then,  even  if  the 
characteristics  of  the  SNOI  change,  the  filter  can  adapt  to  a  new  ‘optimum' solution. 
Such  a  filter  is  referred  to  as  a  ‘Time  Independent  .Adaptive Filter'  (TI.AF);  ‘adaptive’ 
because  the  impulse  response  changes  according  to  the  characteristics  of  the  input, 
and  ‘time  independent’  because  the  changes  in  the  impulse  response  are  not  an 
explicit  function  of  time,  but  rather  a  function  only  of  the  input  to  the  filter.  .An 
example  of  aTlAF  is  the  adaptive  linear  combiner  (ALC)  (10:15-26). 
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The  error  in  a  TIAF  is  a  quadratic  function  of  the  filter  weights;  The  quadratic 
surface  resembles  a  bowl  (10:19,  20).  The  bottom  of  the  bowl  represents  the  .mini¬ 
mum  error,  and  the  filter  strives  to  achieve  that  operating  point.  When  the  statistics 
of  the  input  are  stationary  or  nearly  stationary  (change  very  slowly),  the  filter will  be 
able  to  achieve  and  maintain  a  nearly  optimum  configuration  (^11).  Unfortunately, 
all  digitally  modulated  waveforms  of  interest  in  modern  communication  theory  are 
not  stationary.  They  exhibit  some  periodicity  in  their  statistics  referred  to  as  cyclo- 
stationarity  (2:16-18).  The  TIAF  is  in  general  not  able  to  adapt  quickly  enough  to 
‘track’  the  optimum  solution  that  exists  for  a  cyclostationary  SOI. 

What  is  needed  is  some  way  to  make  the  input  to  the  filter  stationary  so 
that  it  may  better  track  the  optimum  solution.  One  way  to  do  this  when  filtering 
samples  from  a  digitally  modulated  signal  is  to  implement  the  filter  by  using  multiple 
TlAFs  arranged  in  parallel,  and  commutating  the  output  from  each  of  the  TIAFs. 
Then,  each  TIAF  has  its  own  optimum  solution,  or  ‘bowl’,  and  the  filter  can  acheive 
significantly  reduced  overall  MSE  (1:681).  The  selection  of  the  proper  TIAF  is  an 
explicit  function  of  time,  and  hence,  an  adaptive  filter  so  constructed  is  called  a  ‘time 
dependent  adaptive  filter’  (TDAF). 

1.2  Problem  Staiement 

This  thesis  presents  a  TDAF  which  can  be  used  to  improve  the  signal  to  in¬ 
terference  ratio  (SIR)  and  signal  to  noise  ratio  (SNR)  of  digitally  modulated  com¬ 
munications  signals.  The  performance  improvement  of  the  TDAF  over  the  TIAF 
is  determined  based  on  the  application  of  various  metrics,  including  Mean  Square 
Error  (MSE)  and  Bit  Error  Rate  (BER). 

1.3  Scope 

The  intent  of  this  thesis  is  to  produce  a  computer  simulation  of  a  simple  digital 
communications  system  that  can  be  used  to  evaluate  the  performance  of  a  TDAF 
under  realistic  conditions.  A  binary  phase  shift  keyed  (BPSK)  signal  is  used  as  both 
the  SOI  and  the  SNOI.  Comparisons  of  performance  are  made  between  the  TDAF 
and  the  TIAF  for  varying  SNR,  SIR,  interference  carrier  frequency,  and  interference 
data  rate. 
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1.4  Approach 

The  approach  used  to  complete  this  research  is  divided  into  three  phases.  The 
first  phase  has  two  parts.  The  first  part  is  basic  research  into  adaptive  filters  and 
cyclostationary  signals.  The  second  part  is  the  development  of  the  framework  for 
the  computer  simulation  of  the  communication  system.  The  second  phase  is  the 
development  of  the  TDAF  itself,  and  its  insertion  into  the  simulation.  The  last 
phase  is  the  characterization  of  the  TDAF.  Measurement  of  the  TDAFs  ability  to 
accurately  recover  the  signal  is  made  by  running  simulations  for  various  interference 
characterizations. 

1.5  Organization 

Chapter  II  provides  a  discussion  of  the  fundamental  concepts  of  adaptive  filters 
and  cyclostationarity.  Chapter  III  contains  details  on  the  construction  of  the  TDAF 
simulation.  Chapter  IV  presents  the  results  and  analysis  of  data  generated  by  the 
environment.  The  fifth  and  final  chapter  contains  specific  conclusions  along  with 
recommendations  for  future  research. 
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II.  Background 


2.1  Introduction 

In  order  to  understand  how  Time  Dependent  Adaptive  Filters  (TDAF)  im¬ 
prove  the  Signal  to  Noise  Ratio  (SNR)  and  Signal  to  Interference  Ratio  (SIR)  of 
digitally  modulated  bandpass  signals,  there  are  several  concepts  that  need  to  be  un¬ 
derstood.  Among  these  are  the  concepts  of  cyclostationarity  and  adaptive  filtering. 
In  presenting  the  idea  of  cyclostationaritj',  this  chapter  will  provide  the  definition 
of  cyclostationarity,  and  give  an  e.xample  of  a  cyclostationary  signal.  Next,  a  brief 
introduction  to  adaptive  filtering  will  be  presented.  The  Adaptive  Lineai  Combinei 
(ALC)  which  is  a  type  of  Time  Independent  Adaptive  Filter  (TIAF)  will  be  pre¬ 
sented.  The  TIAF  is  presented  first  because  the  TDAF  can  be  implemented  by  com¬ 
bining  multiple  TIAFs  in  parallel.  Finally,  the  concept  of  Time  Dependent  Adaptive 
Filtering  will  be  briefly  introduced.  A  thorough  description  of  the  implementation 
of  the  TDAF  is  reserved  for  Chapter  III. 

2.2  Cyclostationary  Signals 

A  process  is  cyclostationary  if  it  has  a  periodic  components  in  its  autocorre¬ 
lation  function  (2:20).  The  Fourier  series  expansion  of  a  periodic  (with  period  To) 
function  a;(f)  is 

CO 

4t)=  E  (2-1) 

oo 

where 


n 


where  the  coefficient  at  any  given  frequency  is 

C‘n  =  ^  /  (2.2) 

=  (.r(i)e-^’""“”‘)  (2.3) 

Recall  that  the  power  spectral  density  (PSD)  of  any  function  is  the  Fourier  transform 
of  its  autocorrelation.  Therefore,  if  it  can  be  shown  that  the  PSD  of  a  function 
exhibits  delta  functions,  then  the  function  is  cyclostationary. 
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For  ergodic  functions,  the  classical  autocorrelation  function  is  given  by 


Rxx{r)  -  {x{i)x{t  +  t))  (2.4) 

where 


x{t)  =  the  function  of  interest 
T  =  the  lag  value 


Notice  that  it  is  mathematically  equivalent  to  write  Eq  2.4  as 

Rxx{r)  =  {xit  +  T/2)x{t  -  r/2))  (2.5) 

even  though  Eq  2.5  is  not  physically  realizable  (because  of  its  noncausal  nature). 
The  cyclic  correlation  function  is  a  straightforward  extension  of  this  definition 


(2:19) 


nur)  =  {x{l)x(i  +  7)6-^""") 


(2.6) 


where 


x{t)  =  the  function  of  interest 
r  =  the  lag  value 
a  =  the  cycle  frequency 

There  are  two  equivalent  interpretations  of  Eq  2.6.  First,  the  cyclic  autocorrelation 
function  is  the  standard  autocorrelation  but  with  a  time  and  frequency  shifted  ver¬ 
sion  of  itself.  Referring  to  Eq  2.2,  it  can  be  seen  that  a  second  interpretation  is  that 
the  cyclic  autocorrelation  is  the  Fourier  coefficient  of  the  standard  autocorrelation 
at  a  given  cycle  frequency,  a  (2:17-  20).  Notice  that  when  a  =  0,  Eq  2.6  reduces  to 
the  standard  autocorrelation  function  (2:20). 

Wide  Sense  Stationary  (WSS)  random  processes  can  be  either  cyclostationafy 
or  purely  stationary,  but  not  both.  Purely  stationary  processes  are  those  for  which 
no  a  yf  0  can  be  found  to  satisfy  Rx{t)  #  0  except  perhaps  for  the  degenerate  case 
where  r  =  0  (2:20). 


2-2 


If  Eq  2.6  is  re-expressecl  in  the  form  of  Eq  2^5,  and  is  factored  into 

g-j7ror(i+T/2)  g-i5ra(i-T/2)^  interpretation  becomes  apparent.  The  cyclic 

autocorrelation  function  can  be  written  as  a  conventionabcrosscorrelation  ^2:19) 

=  {Hi  +  -  r/2)e-^™(‘-^/2)]) 

=  {nit  +  T/2)v{t  -  r/2))  (2.7) 

=  RHr)  (2.8) 


where 


u{i)  = 

v{t)  =  .a-(<)e+^'^"' 

Assuming  that  at  least  one  of  u{t)  or  v{i)  is  a  zero  mean  process,  the  conven¬ 
tional  cross- correlation  coefficient  is  (8:124) 


Mr)  = 


[r  UMT 


A  normalization  factor  for  the  cyclic  autocorrelation  follows  from  Eq  2.9  and  is  called 
the  temporal  correlation  coefficient  (2:20) 


’  i?x{0) 


(2.10) 


The  following  example  was  taken  from  (2)  and  is  summarized  here  to  show 
how  to  calculate  a  cyclic  autocorrelation. 

Given  a  real  random  purely  stationary  signal  with  zero  mean.  (i{t)  we  can  write 


(a(t))  =  0 


(2.11) 


We  require  that  the  autocorrelation  of  «(t)  be  nonzero: 


{a(i  +  Tf2)a{l  -  r/2))  ^  0 


(2.12) 


2-3 


Since  we  defined  a{t)  to  be  purely  station'll}  (rather  than  cyclostationary),  we  know  that 


{a{t  +  T/2)a{t  —  =  0  for  all  a  0  (2.13) 

Eq  2.13  guarantees  that 

(a(t)e“'’ =  0  for  all  a  51^  0  (2-14) 

Now  consider  the  amplitude-modulated  sinewave 

x(i)  =  a(t)cos{2Tifot  +  9)  (2.15) 

=  i«(t)[e-’'(-"-f''*+'’>  -b  (2.16) 

Multiplying  a{l)  by  cos(27r  fui  -j-  9)  simply  shifts  the  spectrum  of  a{t)  to  ±/o  (and  reduces 
its  magnitude).  Therefore,  since  Eq  2.14  guarantees  that  a(t)  contains  no  finite-strength 
additive  sinusoidal  components,  it  is  clear  that  x{t)  contains  no  finite-strength  additive 
sinusoidal  components  either.  In  other  words,  there  are  not  delta  functions  in  its  power 
spectral  density  function.  If  we  apply  the  non-linear  transformation  that  is  inside  the  time 
averaging  operator  of  Eq  2.5  to  x{t)  we  get 

yr{l)  =  x{t  -b  r/2).r(t  -  r/2) 

=  a{t  -b  r/2)a(t  -  r/2)i[e^-"'^»^  -b  c"-’-'-'''" 

q,^"(-1*/olT29)  _|_  g-7(-1’f/o<+2S)j  (2.17) 

Ne.xt.  we  apply  the  definition  of  the  cyclic  autocorrelation  function  to  t/tCO- 

=  ie^’2--f<'^(«(/,  +  r/2)«(/.  -  r/2)e-^-='“') 

4.ie-J2-/o^(fl(/.  +  r/2)a{t  -  r/2)e-^-"'") 

-b-e'-''(a(t  -b  r/2)a(t  -  r/2)r;-^'-"(°-=^»)‘) 

-bie"-'-®(a(t  -b  r/2)«(t  -  7/2)e“-'‘"^"+-^'’^‘)  (2.1S) 
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Eq  2.13  guarantees  that  the  first  two  terms  of  Eq  2.18  are  zero.  The  last  two  terms  of 
Eq  2.18  can  only  be  nonzero  if  the  exponeni  inside  the  time  averaging  operator  is  zero. 
Clearly,  that  happens  when  a  =  ±2/  -  0.  .Substituting  a  =  ±2/  -  0  into  Eq  2.18,  and 
noting  that  Ra{T)  =  (a(t  +  Tl2)a{i  -  r/2))  we  have  the  cyclic  autocorrelation  function 
function  for  x{t): 


K(r) 


for  a  =  ±2/o 

^  cos(27r /or),  for  a  =  0 

[  0,  otherwise 


(2.19) 


Referring  to  Eq  2.7,  the  correlation  of  the  frequency  shifted  versions  of  .v(i)  represented 
by  u(i)  and  v(i)  should  not  be  surprising.  After  all,  x(i)  was  formed  by  frequency  shifting 
a(i)  by  /o,  and  the  two  spectral  images  of  .r/t)  are  separated  by  2/o. 


Adaptive  Filtering 

2.3.1  Some  Preliminaries.  Given  a  real  V\  SS  random  process  X{t),  we  know 
that  the  autocorrelation  of  X{i)  can  be  written  as  (8:143) 

/?.YA'(r)  =  S{X(0A'(i  +  -)}  (2.20) 


where 


E{-]  =  the  statistical  expectat'K.r  Ofperalor 
T  =  the  lag  value 

If  we  assume  X{t)  is  an  ergodic  process,  then  we  can  restrict  our  attention  to  a  single 
member  function  of  X{t),  x{t)  (8'i7S).  Furthermore,  we  can  sample  x{t)  starting 
at  some  arbitrary  point  k  so  tJiat  we  are  left  with  a  discrete  sequence  .rffc].  For  N 


samples,  this  sequence  can  be  represented  as  a  colurnn  vector  X/.- 


X,  = 


(2.21) 


L  ii-.-v-i) 


forming  an  N  element  array.  Now  the  vale.;  fi  the  a.itororrelation  at  a  specific  lag 
value  can  be  calculated  as  (11) 


iV-l 


k=0 


(2.2?) 


or  equivalently  in  vector  notation 


4’AA  = 


(2.23) 


where 


n  —  the  lag  value 

N  =  the  number  of  samples  in  the  data  vector 
k  =  the  seqence  number  or  vector  element  number 

and  the  superscript  T  represents  the  vc.':lo:-  transpose. 

VVe  can  now  define  the  autocorrelation  matrix  of  a  sampled  sequence 


^*a:[0] 

4>xx[-l]  ••• 

9xx[-n] 

R  = 

<i>xx[A 

^ix[0] 

T  l] 

(2.24) 

4*xx 

^ix{^  1] 

^rx(0] 
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Now  let  us  assuire  tliat  we  are  designing  a  digital  finite  impulse  response  (FIR) 
filter.  Given  the  input  data  vector  we  decide  we  want  the  filter  output  to  be  the 
sequence  </[/c],  called  the  desired  response.  Then  we  can  define  a  ne.'  vector  P,  the 
cross  correlation  between  the  desired  response  and  the  input  vecur  (10:20) 


p  =  E{d\:]Xk}  = 


d[k]Xk 

dlk]Xk-, 

d[k]Xk-iN-i) 


(2.25) 


It  is  now  possible  to  apply  th'-  Wiener-Hopf  equation  to  obtain  the  optimum 
weight  vector  (called  the  Wiener  weight  vector)  (10:22) 


W*  =  R-‘P 


(2.26) 


The  output  of  an  FIR  digital  filter  rjk  is  just  the  convolution  of  the  input  with  the 
filter  weights.  This  can  be  express^'d  in  sampled  sequence  notation  (5:21) 


7/[/o]  =  ^  w(y}’|.t[^’  —  7?] 


(2.27) 


or  equivalently  in  matrix  notation  (10:17) 


r,.  =  W[X, 


(2.28) 


Unfortunately,  neither  R  nor  P  will  in  general  be  known.  In  an  actual  com¬ 
munication  system.,  .c-{k]  will  be  given  by 


.TfA:]  =  d[k]  -f  7(/cj  -f  7i[/;| 


(2.29) 


2-7 


where 


d[k]  =  the  signal  of  interest 
i[/:]  =  some  interfering  signal 
n[k]  =  nois-' 


Substituting  into  Eq  2.22 

=  4  {^^[^'1  +  '^[^i  +  -  ”]  +  -  ”1  +  «(^’  -  ’i]}  (2.30) 

*'  k=0 

In  the  case  of  a  special  sequeicc  r.^  y  be  transmitted  and  anticipated  bj'  the 
receivei,  and  thereby  be  known  a  j  .'ioriai  the  receiver  (7:102),  but  f[^:]  and  n[^]  in 
any  practical  system  are  not  known.  Therefore  it  is  not  possible  to  directly  calculate 
W'. 


S.3.2  The  Least  Mean  Squared  (LMS)  AUjoriihm.  If  the  input  to  the  filter  is 
stationary,  and  W  is  allowed  to  vary  according  to  some  rule  which  tends  to  mini¬ 
mize  the  difference  between  the  actual  filter  response  and  the  desired  response,  then 
W  becomes  an  arbitrarily  close  approximation  of  W'  (11).  The  equation  for  the 
difference  between  d[k]  and  the  filter  response  y\k]  is  simply 

e|(-|  =  rfiA-l  -  jl/n]  (2.31) 

Note  that  the  power  in  £(/.:]  is  the  Mean  Squared  Enor  (MSE)  and  is  given  by 

f  =  (2-32) 

It  can  be  shown  (10:19,20)  that  ^  is  a  quadratic  function  of  the  filter  tap 
weights,  the  desired  res  joiise,  and  the  filler  input.  ^  is  referred  to  as  the  performance 
surface,  and  for  a  two  tap  filter,  it  is  a  parabloid  (a  Ik  perparabloid  if  there  are  more 
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Figure  2.1.  A  quadratic  mean  squared  error  surface.  The  verticle  axis  is  (  and  the 
horizontal  axes  are  the  filter  weights 

than  2  weights)  and  resembles  a  bowl  (see  Figure  2.1).  The  bottom  of  the  bowl  is 
the  minimum  MSE  and  when  the  bowl  is  projected  onto  the  weight-vector  plane, 
the  minimum  describes  the  point  in  weight  space  described  by  W’  (10:21). 

Given  that  d[k]  is  Imown,  it  is  possible  to  search  the  error  surface  for  its  mini¬ 
mum,  thereby  arriving  at  the  optimum  weight  vector  W”.  The  most  straightforward 
method  of  searching  the  error  surface  is  the  LMS  algorithm  (7:1 01 ).  Other  algorithms 
exist  and  are  presented  in  (10:Chap  S). 

A.  complete  derivation  of  the  LMS  algorithm  is  given  in  Adaptive  Signal  Pro¬ 
cessing  by  VV’idrow  and  Stearns  (10:99-101);  only  the  result  is  given  here: 


Wa.+,  =W,.-f2/te,.X, 


(2..33) 


where 


//  =  the  adaptation  coefficient 


In  words,  the  weight  vector  for  the  next  input  vector  is  the  sum  of  the  current  weight 
vector  and  a  scaled  product  of  the  error  and  the  current  input  vector. 

Eqs  2.27,  2.31,  and  2.33  define  an  adaptive  filter  incorporating  the  LMS  algo¬ 
rithm.  The  two  important  cissumptions  here  are  first,  that  an  acceptable  estimate  of 
the  desired  signal  dk  is  available,  and  second  that  Xa-  is  a  stationary  process.  Should 
Xfc  fail  to  be  a  stationary  process,  then  the  bottom  of  the  bowl  defined  by  the  the 
surface  ^  will  exist  at  different  points  in  weight  space.  The  statistics  of  Xa-  must 
be  non-varying  or  vary  slowly,  so  that  W  will  in  general  approach  a  near  optimum 
solution  (11). 

Extending  the  bowl  analogy,  searching  ^  for  its  minimum  is  equivalent  to  rolling 
a  marble  down  the  side  of  the  bowl.  The  steepness  of  the  sides  of  the  bowl  is 
determined  by  the  power  in  d[k]  relative  to  the  power  in  n[k]  and  ?'[/;]  (11).  The 
more  power  in  d{k],  the  steeper  the  sides  of  the  bowl;  hence,  the  more  rapid  the 
convergence. 


With  adaptive  filters,  given  an  error  surface  there  is  a  trade-off  between  rapid 
convergence  and  close  approximation  of  W*.  The  trade-off  involves  selection  of  an 
appropriate  adaptation  coefficient,  //.  The  smaller  /q  the  more  closely  W  approaches 
W*  (assuming  stationarity  of  the  input).  For  faster  adaptation,  /t  is  chosen  to  be 
larger,  but  the  MSE  also  increases.  A  good  rule  of  thumb  is  to  select  such  that 
(10:103,111-  114) 


''  {N  -}-  l)(Power  in  Xt)  ^ 

where  M  is  the  misadjustmevL  of  the  filt.,r.  The  misadjustment  is  a  measure  of  the 
average  distance  between  W  and  W’  that  the  filter  designer  is  willing  to  live  with. 
Smaller  M  results  in  close  approximation  of  W',  but  slower  adaptation.  Eq  2.34 
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also  assures  convergence  of  the  filter  for  0  <  M  <  1.  Unlike  nonadaptive  BIRs,  the 
adaptive  transversal  filter  can  fail  to  converge  (10:102). 


2.3.3  Time  Dependent  Adaptive  Filteiing.  Digitally  modulated  signals  are 
not  stationary,  but  rather  cyclostationary  (4:1).  As  a  result,  there  does  not  exist  a 
single  value  of  W*  associated  with  an  single  error  surface  Instead,  there  exist  (for 
a  sampled  signal)  a  finite  number  of  error  surfaces  (1:679,  680).  A  TDAF  simply 
provides  a  separate  TIAF  for  each  error  surface  in  weight  space.  This  can  result  in 
a  significant  reduction  in  the  MSE  of  the  filter  output  (7:3).  Each  TIAF  in  the  TDAF 
has  an  independently  adapted  weight  vector,  W^.^.  The  TDAF  LMS  algorithm  is 
(1:681) 

Vk  =  (2.35) 


and 


W 


+  2pCkXk  Ck  due  to 

< 

otherwise 


(2.36) 


For  error  surface  the  weight  vector  W/.j  is  convolved  with  the  input  and 
is  subsequently  updated  according  to  the  LMS  algorithm.  The  other  weight  vectors 
Wkjfy  are  dormant.  Hence,  the  overall  time  to  adaptation  for  the  TDAF  is  slower 
than  for  the  TI.AF  since  each  weight  vector  is  updated  only  periodically,  rather  than 
at  each  sample  time  k  as  in  the  case  of  the  TI.AF.  The  specific  implementation  of 
the  TD.AF  will  be  covered  in  detail  in  Chapter  III. 


2..^  Chapter  Summai-y 

This  chapter  provided  an  overview  of  the  concepts  of  cyclostationarity  and 
adaptive  filtering.  A  mathematical  means  by  which  to  determine  if  a  signal  is  cy- 
clostationary  was  provided,  and  an  example  of  the  calculation  was  presented.  Time 
Independent  Adaptive  Filtering  was  covered  to  form  a  foundation  for  the  introduc¬ 
tion  of  the  concept  of  Time  Dependent  Adaptive  Filtering. 
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III.  Simulation  Implementation  and  Verification 

3. 1  Introduction 

This  chapter  provides  a  detailed  description  of  the  methodology  employed  to 
characterize  the  performance  of  the  Time  Dependent  Adaptive  Filter  designed  in 
support  of  this  research  effort.  Included  is  a  description  of  the  computer  program 
that  was  designed  and  written  to  facilitate  the  characterization. 

3.2  Top  Level  Description  oI  the  Simulation 

The  simulation  is  a  computer  program  that  was  written  in  the  C  programming 
language,  and  compliled  on  the  public-domain  GI^’U  C  compiler.  Pains  were  taken 
to  adhere  to  the  .ANSI  standard  that  was  recently  established  for  the  C  language. 
.\s  a  result,  the  code  should  compile  and  run  on  any  system  that  has  an  ANSI  C 
compiler  and  sufficient  memory  resources,  including  an  IBM  AT  class  compatible 
computer.  In  fact,  much  of  the  code  was  developed  on  a  PC  compatible  computer 
with  an  S03S6  microprocessor  and  80-387  floating  point  coprocessor  using  Borland 
C-r+  in  .ANSI  mode. 

The  program  simulates  a  simple  communication  system.  Included  in  the  system 
are  a  data  source,  bandlimiting  filters,  a  modulator,  a  noisy  channel,  adaptive  filters 
and  a  demodulator  (See  Figure  3.1).  The  SOI  cjiannel  and  SNOI  channel  shown  in 
Figure  3.1  are  identical.  The  data  rate  for  the  SNOI  channel  can  be  set  to  some 
percentage  of  the  SOI  data  rate. 

In  the  simulation,  the  .symbol  rale  of  the  signal  of  interest  (SOI)  channel  is 
1  symbol/second.  All  other  times,  rates,  and  frequencies  are  based  on  that  value. 
While  an  actual  practical  communication  system  is  unlikely  to  use  such  a  data  rale, 
it  is  a  simple  matter  to  scale  that  rate  to  the  appropriate  level.  This  normalization  to 
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Figure  3.1.  Signal  flow  through  the  simulation 
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the  SOI  data  rate  simplifies  the  design  of  the  simulation  and  analysis  of  the  resulting 
data. 

There  is  a  fundamental  “unit  of  simulation  time”  referred  to  as  an  epoch.  An 
epoch  always  is  comprised  of  an  integer  number  of  symbols.  The  duration  of  the 
simulation  is  determined  by  specifying  the  number  of  symbols  per  epoch  and  the 
number  of  epochs  in  the  simulation. 

Every  reasonable  attempt  has  been  made  to  modularize  the  software  to  fa¬ 
cilitate  modification  of  the  simulation.  Each  component  of  the  simulation  (filters, 
modulators,  etc.)  accepts  a  single  value  for  the  input  and  provides  back  to  the  pro¬ 
gram  a  single  value  as  output.  Hence,  it  is  straightforward  to  substitute  a  different 
filter  or  modulator  for  the  one  provided. 

3.2.1  BER  Version  and  LC  Version  There  are  two  distinct  versions  of  the 
simulation  which  differ  in  the  data  that  they  produce,  but  not  in  the  filtering  of  the 
composite  channel  signal.  The  first  version  allows  the  filters  to  adapt  for  a  specified 
number  of  epochs,  and  then  freezes  the  weights  and  continues  to  run,  keeping  track 
of  the  number  of  bit  errors  for  the  TDAF  channel  and  for  the  TIAF  channel.  This 
is  referred  to  as  the  BER  (Bit  Error  Rate)  Version. 

The  second  version  allows  the  filters  to  adapt  continuously  while  storing  the 
resulting  squared  error.  .After  a  full  epoch,  the  filler  weights  are  reinitialized  and 
again  allowed  to  adapt.  The  resulting  MSE  .  averaged  with  that  of  the  previous 
epoch.  The  process  is  repeated  a  specified  number  of  epochs  (selectable  at  run  time) 
providing  an  estimate  of  the  expected  value  of  the  MSE  as  a  function  of  adaptation 
time.  This  is  called  the  LC  (Learning  Curve)  Version. 

3.2.2  Input,  Output,  Notation,  and  Definitions.  Depending  on  the  version 
being  run,  the  simulation  requires  24  input  parameters  at  run  time.  The  required 
parameters  can  be  ,^ed  in  at  the  terminal,  but  ease  of  use  is  significantly  enhanced 
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if  the  parameters  are  typed  into  a  text  file,  and  provided  to  the  simulation  using  the 
UNIX  or  DOS  redirection  operator.  The  parameters  required  for  the  BER  version 
are  listed  in  ijVppendix  A  along  with  a  list  of  the  data  files  produced.  The  same 
information  for  the  LC  version  is  listed  in  Appendix  B. 

The  MSE  returned  by  the  BER  version  of  the  program  is  different  from  the 
MSE  returned  by  the  LC  version.  The  BER  version  of  the  program  freezes  the  filter 
weights,  and  ^hen  starts  to  collect  squared  error  data.  That  data  is  then  averaged 
and  returned  by  the  program  as  the  MSE.  The  LC  version  allows  the  filter  weights 
to  continuously  adapt,  and  averages  squared  error  sample  by  sample  for  all  epochs 
yielding  a  vector  of  MSE  .is  a  function  of  sample  number.  The  MSE  error  for  the 
last  P  symbols  of  the  learning  curve  is  averaged,  and  that  number  is  returned  by 
the  program  as  the  MSE.  As  a  result,  for  equal  inputs  and  initialization,  the  MSE 
returned  by  the  LC  version  of  the  program  will  always  be  smaller  than  that  returned 
by  the  BER  version. 

All  points  in  the  signal  path  of  the  simulation  are  necessarily  discrete  samples. 
In  an  actual  system,  many  of  these  signal  would  be  continuous  time.  In  this  chapter, 
continuous  time  signals  are  represented  as  functions  of  time:  x{t),  etc.  Discrete 
time  signals  are  represented  as  functions  of  sample  number:  or  Xk,  f/(A']  or  (k-^ 

etc.  Here  is  a  list  of  most  of  the  signals  used  in  the  simulation: 

•  ^S/VO/(0-  baseband  data  .signal  (1  or  0). 

•  wso/(0’  ^S/VO/(0'  formatted  baseband  data  signal  (bi-polar  or  bi-phase). 

•  bandlimited  ba,seband  data  signal. 

•  dsoi{t),  dsNoi{i)-  modulated  signal. 

•  n{t):  Additive  white  gaussian  noise. 

•  x{t):  sum  of  dsoiit),  dsNoi{i) 
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•  yfDAFii)^  yTiAF{t)-  adaptively  filtered  bandpass  signal. 

•  y'rDAFi^)^  y'TiAFi^)-  demodulated  signal. 

•  '>^TDAF{t),  recovered  signal. 

•  ^TDAF{i)i  ^TiAF{t)'-  difference  signal,  d{t)  —  y{t) 


3.3  Simulation  Details 

The  simulation  has  two  time  bases:  one  for  the  SOI  and  another  for  the  SNOI. 
The  simulation  operates  by  calculating  the  signal  value  at  each  point  in  the  signal 
path  and  then  incrementing  time  by  l/(sample  frequency).  The  signal  value  is  then 
recalculated.  This  process  is  repeated  until  the  required  number  of  data  symbols 
have  been  processed. 

3.3.1  Data  Generator.  The  data  generator  (DG)  used  in  this  simulation  pro¬ 
duces  a  pseudorandom  sequence  of  I’s  and  O’s  sampled  at  the  appropriate  frequency. 
Incorporated  in  the  DG  is  a  routine  taken  directly  from  Numerical  Recipes  in  C 
(6:220,  228)  called  irbit2  which  actually  calculates  the  value  of  each  data. bit.  When¬ 
ever  a  full  bit  time  has  elapsed,  the  DG  calls  irbit2,  requesting  a  new  bit  value,  either 
1  or  0,  for  the  random  sequence.  The  DG  continues  to  output  that  value  for  the 
next  full  bit  time. 

The  routine  irbit2  lequires  a  primitive  polynomial  modulo  2,  referred  to  here 
as  the  seed,  to  form  the  random  bits.  The  polynomial  used  was 

.t’®  +  +  x^  +  .T°.  (3.1) 

In  C,  the  above  polynomial  is  represented  as  the  bit  sequence 


{lOOOOOOOOOOOOiOOlll} 


(3.2) 


and  generates  a  random  sequence  with  a  period  of  repetition  of  2^®  — 1  =  262143  bits. 
The  longest  simulation  run  was  512  epochs  x  16  symbols  per  epoch  =  8192  symbols. 
Hence,  as  far  as  the  simulation  was  concerned,  the  data  sequence  was  purely  random. 

A  call  to  the  DG  routine  requires  six  input  parameters: 

•  time:  The  current  value  of  the  base  clock  for  the  signal  path  under  considera¬ 
tion  (either  SOI  or  SNOI). 

•  datarate:  This  is  constrained  to  be  1  for  the  SOI.  For  the  SNOI,  it  is  calcu¬ 
lated  by  the  simulation  as  1/SNOI  symbol  frequency  where  the  SNOI  symbol 
frequency  is  an  input  parameter  to  the  program. 

•  lastdata:  The  value  of  the  last  sample  returned  by  DG.  This  variable  is  modified 
by  the  routine. 

•  lasttime:  The  time  when  the  current  bit  started. 

•  iseed:  The  current  value  of  the  .seed  required  by  irbitS.  This  variable  is -modified 
by  the  routine. 

•  oulpxUflag:  A  flag  which  forces  the  DG  to  return  a  square  wave  rather  than 
random  data  when  it  is  set. 


3.S.S  Finite  Impulse  Response  Loxo  Pass  Filters.  FIR  filters  are  used  for  the 
required  low  pass  filters  (LPF)  because  the  filter  coefficients  can  be  quickly  and 
easily  calculated  so  that  it  is  possible  to  select  the  cutoff  frequency  of  the  filter 
at  run  time.  The  FIR  filter  coefficients  are  calculated  using  the  windowed  Fourier 
Transform  method  with  the  Hamming  window(5;444-452): 


hk  = 


XUk  = 


sin[u;c(^:  —  •/V/2)] 
ir{k  -  N/2) 

I  0.54  —  0.46  cos(27r /:///),  0  <  k  <  N 
I  0,  otherwise 


(3.3) 

(3.4) 
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where 


Wo  =  the  discrete  cutoff  frequency  (0  <  <  tt) 

N  =  the  number  of  taps  in  the  FIR  filter 

The  output  of  the  filtering  routine  is  the  convolution  of  the  filter  coefficents  given 
by  Eq  (3.3)  and  the  past  N  input  values: 

/V-l 

ni'[k]  =  ^  h[n]m[k  —  n]  (3.5) 

n=0 

The  filters  are  used  for  two  purposes  at  two  different  points.  First,  to  band 
limit  the  baseband  message  signal  (either  SOI  or  SNOI)  at  the  input  to  the  channel, 
and  second,  to  reject  the  double  frequency  component  of  the  demodulated  signal  at 
the  output  of  the  mixer  in  the  demodulator  (see  Section  3.3.4). 

A  call  to  the  LPF  routine  requires  five  input  parameters: 

•  input:  The  value  of  the  signal  to  be  filtered. 

•  niimdelays:  One  less  than  the  number  of  filter  coefficients. 

•  weights:  A  vector  containing  the  values  of  the  filter  coefficients. 

•  contents:  A  vector  containing  the  past  numdelays  samples  of  the  input. 

•  gain:  The  output  of  the  filter  is  simply  multiplied  by  this  number. 

3.3.3  Data  Formatter  The  data  formatter  (DF)  simply  converts  the  binary 
digits  of  the  data  generator  into  a  bi  polar  signal.  When  its  input  is  a  1,  the  formatter 
does  nothing.  When  its  input  is  a  0,  the  formatter  returns  -1.  The,  formatter  is  also 
capable  of  formatting  the  output  into  a  Manchester  code.  However,  no  data  was 
taken  with  the  formatter  running  in  this  mode. 
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3.3.4  Modulator  and  Demodulator  BPSK  modulators  have  several  advan¬ 
tages  that  lend  themselves  to  use  in  a  simulation  such  as  this.  First,  they  are  simple 
to  code.  Second,  they  are  in  common  use.  Finally,  the  demodulated  signal,  once  low 
pass  filtered,  is  restored  to  the  shape  of  the  baseband  signal.  As  a  result,  a  BPSK 
modulator  was  chosen  for  this  simulation. 

The  equation  that  describes  the  modulation  process  is  (9:130) 

5(Z)  =  Am'{t)  cos{2ir  Jet)  (3.6) 

where 

A  =  carrier  amplitude 
fc  =  carrier  frequency 
m'{i)  =  the  formatted  message 

Note  that  if  m'(t)  =  ±1,  the  power  in  the  modulated  signal  is  (9:16) 

P=Y  (3'^) 

A  call  to  the  modulator  routine  requires  four  input  parameters: 

•  input:  The  value  of  the  baseband  signal. 

•  carrierfreq:  The  carrier  frequency,  fc- 

•  carricrampl:  The  carrier  amplitude,  A. 

•  time:  The  current  value  of  the  base  clock  for  the  signal  path  undei  considera¬ 
tion. 
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The  demodulator  simply  restores  the  bandpass  signal  to  baseband.  The  equa¬ 
tion  that  describes  its  operation  is 

ih{i)  =  y{t)  cos(2^ fci  -b  o)  (3.8) 

The  output  of  the  demodulator  is  fed  to  a  low  pass  filter  to  reject  the  double  frequency 
component  that  results  from  Eq  3.8. 

A  call  to  the  demodulator  routine  requires  four  input  parameters: 

•  input:  The  value  of  the  bandpass  signal. 

•  carrierfreq:  The  carrier  frequency,  fc- 

•  phase:  A  correction  factor  provided  in  the  event  that  there  is  a  significant 
phase  shift  through  the  channel. 

•  time:  The  current  value  of  the  base  clock  for  the  signal  path  under  considera¬ 
tion. 

3.3.0  Noise  Generator  The  signal  n(t)  is  provided  by  a  routine  called  gasdev 
taken  from  Numerical  Recipes  in  C  (6:210.  211,  217).  The  routine  gasdev  returns 
normally  (Gaussian)  distributed  random  samples  with  zero  mean  and  unit  variance. 
The  simulation  multiplies  the  returned  sample  b\  a  gain  factor  in  oider  achieve  noise 
power  other  than  0  dBW.  As  an  example,  assume  the  desired  level  of  noise  jjower  is 
-1-5  dBW.  Converting  to  absolute  power  levels 

-  1.5  dBW  =  iQ-’®/’*’  =  0.03162  watts  (3.9) 


A  normally  distributed  random  process  X  with  mean  n  and  variance  can  be 
converted  into  a  normal  random  process  Z  with  zero  mean  and  unit  variance  by 
applying 


X-ti 


<7 


(3.10) 


3-9 


Therefore,  Z  can  be  converted  into  A'"  by  applying 

X  =  aZ^  /t  (3.ir) 

Since  is  the  ac  power  in  the  process  (we  want  to  maintain  zero  mean,  hence  zero 
dc),  multiplying  Z  by  the  square  root  of  the  desired  power  yields  the  appropriate 
random  process.  So  in  the  example, 

,7^  =  0.03162 
a  =  0.17783 

Providing  the  simulation  with  a  noise  gain  of  0.17733  at  run  time  results  in  -15  dBVV 
of  noise  ])ower  being  added  to  the  SOI  and  SNOJ. 

The  only  input  parameter  gasdev  requires  from  the  simulation  is  a  seed  for 
the  random  number  generator.  Any  nega.tive  integer  on  the  interval  [-65536,  -1]  is 
acceptable. 

2.3.6  Time  Independent  Adaptive  Filter  (TIAF)  The  TIAF  is  adapted  from 
Chapter  6  of  VVidrow  and  Stearns  Adaptive  Signal  Proeessing  (10).  The  transver¬ 
sal  adaptive  linear  combiner  (ALC)  has  the  advantage  that  it  is  easy  to  code  and 
the  .single  input,  single  output  characteristic  of  the  filter  is  perfectly  suited  to  this 
simulation,  A  block  diagram  of  a  transvensal  adaptive  linear  combiner  is  shown  in 
Figure  3.2.  Note  the  presence  of  a  bias  weight  in  Figure  3.2.  This  allows  for  more 
rapid  convergence  of  the  filter  when  some  dc  or  very  low  frequency  component  is 
present  in  the  signal  being  filtered  (1:679). 

The  TIAF  routine  operates  on  two  vectors  which  are  passed  in  to  the  routine 
each  time  it  is  called.  The  first  is  the  vector  W^.  containing  the  filter  coefficients. 
The  second  is  the  input  vector  X*  that  contains  the  past  N  inputs  to  the  filter, 
where  N  is  the  number  of  filter  coefficients. 
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Figure  3.2.  An  Adaptive  Linear  Combiner 

The  TIAF  is  initialized  by  calling  a  routine  that  sets  Wjt  =  Xk  =  0.  In  the 
HFR  version  of  the  simulation,  the  TIAF  is  initialized. only  once.  In  the  LC  version, 
i.ne  TIAF  is  re-initialized  after  evv.y  P  .symbols,  where  P  is  the  number  of  symbols 
in  one  epoch. 

A  call  to  the  TI.^F  routi  le  requires  eight  input  parameters: 

•  input:  The  signal  being  filtered. 

•  desired:  The  training  signal. 

•  mu:  The  adaptation  coefficient. 

•  numtaps:  The  number  of  filter  coefficients. 

•  error:  The  value  of  dk  —  j//...  Recall  that  dk  is  the  desired  signal  and  ijk  is  the 
value  returned  by  the  routine.  This  variable  is  modified  by  the  routine. 

•  weights:  The  vector  W/..  representing  the  filter  weights.  This  vector  is  modified 
by  the  routine. 


3-11 


•  contents:  The  vector  X/.-  representing  the  past  inputs  to  the  filter.  This  vector 
is  modified  by  the  routine. 

•  bias:  The  value  of  the  bias  filter  weight.  This  variable  is  modified  by  the 
routine. 

When  the  TIAF  routine  is  called,  the  process  illustrated  in  Figure  3.3  is  per¬ 
formed. 


3.3.7  Time  Dependent  Adaptive  Filter  (TDAF)  The  signal  m'[i)  defined  in- 
Section  3.2.2  meets  the  criteria  for  d{t)  defined  by  Eqs  2.11,  2.12  and  2.13  inSec- 
tion  2.2.  Therefore  the  output  of  the  modulator  ofSection  3.3.4,  y{t),  is  a  cyclosta¬ 
tionary  signal  with  cycle  frequencies  db2/c. 

The  TDAF  is  made  up  of  a  number  of  TIAFs  with  a  common  input.  The 
output  from  the  TDAF  is  commutated  from  the  output  of  the  TIAFs  as  shown  in 
Figure  3.4.  The  number  of  TIAFs  used  is  equal  to  the  number  of  samples  taken  per 
symbol  of  the  SOI.  Hence,  a  priori  knowledge  of  the  sampling  rate  and  symbol  rate  of 
the  SOI  is  required  to  implement  this  filter.  Note  that  the  memory  requirements  for 
the  TD.AF  are  larger  than  that  for  the  TIAF,  because  of  the  multiple  weight  vectors 
required.  The  processing  time  is  approximately  the  same  for  both  filters  however, 
since  each  call  to  cither  routine  involves  only  a  single  convolution. 

To  show  how  a  TD.AF  stationarizes  a  cyclostationary  signal,  consider  the  time 
domain  signal  consisting  of  a  periodic  pulse  train  with  period  To  and  sampling  fre¬ 
quency  1/Ta  (see  Figure  3.5): 

L-\ 

x{i)  =  6(1-  nTo)  (3.12) 

71=:— L 


Let 


v{t,T)  =  x{t)x{i  -b  t) 


x{t),  T  =  kTo]  k  =  0,  ±1,  ±2,  ■  ■  • 
0,  otherwise 


(3.13) 
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{Initialize  v<ariable  that  will  hold  output  value... 
sum  =  0.0; 

(An  N  element  vector  ranges  from  Xo  to  x^j^t.j^... 

(Shift  the  input  vector  X|j  to  the  right  and  multiply  input  by 
filter  weights  while  summing  the  product... 

for  count  =  (N  -  1)  to  1  do 
begin 

^COUIlt  “  ^COllIlt— I: 

sum  =  sum  +  x  count  *  count; 
end; 

(now  take  care  of  the  current  input 

Xo  =  input; 

sum  =  sum  +  xo  *  wq; 

(now  take  care  of  the  bias  input 

sum  =  sum  +  1  *  wijja,; 


(The  current  value  of  sum  will  be  returned  to  the  program  as  the 
output... 

(but  the  error  must  still  be  calculated 
error  =  desired  -  sum; 

{Apply  the  LMS  algorithm  to  the  filter  weights... 

for  count  =  0  to  N  —  1  do 
begin 

^*'count  —  ^^count  ■^  ‘2  *  ji  *  error  *  .'Tcounli 
end; 

return  sum; 


Figure  3.3.  Pseudocode  for  the  TIAF  algorithm 


Figure  3.4.  A  TD.AF  implemented  as  a  parallel  bank  of  TIAFs 
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Figure  3.5.  The  impulse  train  described  by  Eq  3J2  sampled  every 


From  Eq  2.6  we  have 

/?(!■)  =  (3.1.1) 

=  I'm  ^  [  j(i,  (3.15) 

To— oo  Jq  JTo 

—  hm  ^1  y{Lr)[cosi2nat)  —  jsin{2t:cti.)](ll  (3.16) 

7b— •»  Jq  JTo 

—  Jr  /  y(i-r)cos(2»al)(Il  —  ^  f  y{L.T)s\n{27ToL)dl  (3.17) 

ra-«>  Tq  JTo  '  To  JTo 

Recall  (Section  2.2)  that  a  signal  is  cyclostationary  if  J^(t)  ^  0  for  some  a  ^  0. 
Substituting  Eq  3.12  into  Eq  3.17  and  picking  a  =  l/J’o.  ’fter  interchanging  the 
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order  of  summation  and  integration,  we  have 

=  i  f;  f Ml- n%)c.os{iM)di- 
?’  r  1 

TfrY,  S{t—,nTo)sm{2Ti—i)di  (3.18) 

^On=-L'''^'o  -^0 

=  ^  E  W27r;^n'io) -isin(27rinro)  (3.19) 

■^0„=_£  -‘o  iQ 

=  2L/7o  (3.20) 

I 

which  is  clearly  not  equal  to  zero.  Therefore  Eq  3.12  defines  a  cyclostationary  signal. 

Now  consider  that  x(L)  is  applied  to  the  filter  of  Figure  3.4.  In  the  first  TIAF 
of  the  bar'-  the  pulse  occurs  at  the  first  delay  element;  in  the  second  TIAF,  the 
pulse  occurs  at  the  second  delay  element,  and  so  on.  After  a  full  period  of  the  signal, 
the  next  pulse  in  the  train  will  be  applied  to  the  filter,  where  it  will  occur  at  the  first 
delay  element  of  the  first  TIAF.  See  Figure  3.3.7.  Thus,  at  the  time  that  any  given 
TIAF  in  the  TDAF  is  updated,  the  pulse  will  always  be  in  the  same  bin.  Therefore, 
from  the  perspective  of  the  Nth  TIAF,  the  signal  can  be  written  as  the  nonperiodic 
function 

^TIAF{t)  -  h{t  -  N)  (3.21) 

The  autocorrelation  of  .t’TMF(0  exists  only  foi  r  =  0.  Therefore,  except  for  the. 

degenerate  case  of  r  =  0  the  cyclic  autocorrelation  is  zero.  Hence,  the  TIAF  sees  a 

purely  stationary  signal. 

.*\  call  to  the  TDAF  routine  requires  nine  input  parameters: 

•  mput:  The  signal  being  filtered. 

•  desired:  The  training  signal. 

•  mu:  The  adaptation  coefficient. 
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TIAF  0  at  t  =  io 
TIAF  \  Sit  t  =  to  Ts 
TIAF  2  at  t  =  ^0  +  2Tg 
TIAF  3  at  i  =  to  + 

TIAF  4  at  t  =  to  "F 

TIAF  0  at  t  =  to  4"  =  to  +  2o 

TIAF  latt  =  to+T^+2o 


Figure  3.6.  Stationarization  of  a  periodic  pulse  train  by  a  TDAF 
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•  sampersym:  The  number  of  samples  per  symbol  of  m'{t)  (also  the  number  of 
TIAFs  in  the  TDAF). 

•  nunitaps:  The  number  of  filter  coefficients  in  each  TIAF. 

•  tdafwts:  The  two  dimensional  array  representing  the  filter  weights.  This 
array  is  modified  by  the  routine. 

•  contents:  The  vector  Xfc  representing  the  past  inputs  to  the  filter.  This  vector 
is  modified  by  the  routine. 

•  error:  The  array  containing  the  values  of  dk  —  pk-  Each  TIAF  has  an  associated 
error,  hence,  the  error  must  be  stored  in  an  array.  This  array  is  modified  by 
the  routine. 

•  bias:  The  array  containing  the  values  of  the  bias  filter  weights.  Each  TIAF 
has  an  associated  bias  weight.  This  variable  is  modified  by  the  routine. 

Whenever  the  TDAF  routine  is  called,  the  process  illustrated  in  Figure  3.7  is 
performed.  Note  the  similarity  to  Figure  3.3. 

3Jf  Code  Verification 

When  properly  coded,  the  adaptive  filters  in  this  simulation  converge  to  a  zero 
mean  squared  solution  if  there  is  no  interference  or  noise  in  the  input  signal  (i.e. 
when  =  d{t))  (11).  Figure  3.8  shows  the  converged  response  of  the  filters  under 
those  circumstances. 

The  small  amount  of  ripple  visible  in  the  response  of  the  TDAF  is  due  to  its 
slow  convergence.  This  was  verified  by  decreasing  the  adaptation  time  and  noting 
increased  ripple.  Each  input  of  the  TIAF  gives  rise  to  an  adjustment  of  the  filter 
weights.  For  the  TDAF,  however,  the  individual  TIAFs  comprising  the  filters  are 
adjusted  only  once  per  symbol.  Therefore,  the  TDAF  converges  much  more  slowly. 


} 


{Initialize  variable  that  will  hold  output  value... 
sum  =  0.0; 

{If  this  is  the  n*'‘  call  to  the  routine,  then  the  TIAF  that  is  active  is...  } 

tiafnum  =  (n  mod  sampersym)  +  1; 

{Increment  n  for  the  next  call  to  the  routine  } 

n  =  n  +  1; 

{Shift  the  input  vector  X|.  to  the  right  and  multiply  input  by 

filter  weights  while  summing  the  product...  } 

for  count  =  (N  -  1)  to  1  do 
begin 

^count  —  ^count-1) 

sum  =  sum  +  .x  count  *  count,  tiafnumi 
end; 

{now  take  care  of  the  current  input  } 

.Xo  =  input; 

sum  =  sum  +  Xq  *  w^^  tiafiiumi 

{now  take  care  of  the  bias  input  } 

sum  =  sum  +  1  *  tiafnumi 

{The  current  value  of  sum  will  be  returned  to  the  program  as  the 

output,  but  the  error  must  still  be  calculated...  } 

tiafnum  =  <^esired  -  sum; 

{...and  the  filter  weights  must  be  updated...  } 

for  count  =  0  to  A  —  1  do 
begin 

’^count,  tiafnum  ~  ^^count,  tiafnum  +  2  *  /x  +  error  tiafnum  *  •'’'count! 
end; 

return  sum; 

I 


Figure  3.7.  Pseudocode  for  the  TDAP  algorithm 


Figure  3.8.  Convergence  of  the  adaptive  filters  with  no  interference  or  noise;  a)  the 
desired  signal;  b)  the  TIAF  response;  c)  the  TDAF  response. 
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3.5  Chapter  Summary 

This  chapter  provided  a  detailed  description  of  the  simulation  developed  to 
compare  the  performance  of  a  Time  Dependent  AdvHptive  Filter  with  that  of  a  Time 
Independent  Adaptive  Filter.  In  the  simulation,  both  the  TDAF  and  TIAF  simul¬ 
taneously  filter  the  channel  signal.  The  simulation  calculates  the  resultant  Mean 
Squared  Error  for  both  filters,  another  version  of  the  simulation  also  counts  the 
number  of  bit  errors  at  the  receiver  end. 
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IV.  Results 


4.1  Introduction 

In  this  chapter,  the  results  obtained  from  a  number  of  simulation  runs  are 
presented.  Each  simulation  run  varied  one  or  more  of  the  following  parameters: 

•  Signal  to  Noise  Ratio  (SNR)  at  the  filter  input 

•  Signal  to  Interference  Ratio  (SIR)  at  the  filter  input 

•  Baud  rate  of  interferer 

•  Carrier  frequency  of  the  interferer 

The  improvement  factor  for  the  TDAF  is  defined  as  (7:51) 

_  MSEtiaf 
MSEjoaf 

The  improvement  factor  will  be  used  as  the  primary  figure  of  merit  for  the  TDAF. 
A  comparison  of  Bit  Error  Rate  (BER)  is  also  made,  but  it  is  emphasized  that 
the  results  presented  are  only  very  rough  estimates  of  the  actual  BER  because  only 
relatively  short  (8192  symbols)  simulations  were  run.  Longer  simulation  runs  were 
po.ssible,  but  with  a  sample  frequency  of  64  samples  pei  symbol,  the  run  time  involved 
would  have  reduced  the  number  of  runs  that  could  be  made.  Furthermore,  computer 
disk  space  was  extremely  scarce,  and  the  fact  that  longei  runs  result  in  longer  data 
files  contributed  to  the  decision  to  limit  the  length  of  each  simulation  run. 

Simulation  in  a  Noisy  Environment 

The  first  runs  of  the  simulations  were  made  with  with  the  carrier  amplitude 
of  the  interfering  signal  set  to  zero.  The  amplitude  of  the  SOI  carrier  and  the  noise 
gain  are  shown  in  Table  4.1.  The  remaining  input  parameters  were  not  changed  from 
one  simulation  to  the  next.  Their  values  are  shown  in  Tables  4.2  and  4.3. 
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Table  4.1.  SOI  carrier  amplitude  and  noise  gain  for  noisy  environment  simulations 


Parameter 


random  bit  generator  seed 


.AVVGN  seed 


number  of  samples/symbol 


Manchester  or  Bipolar  format  (l=Bi,  0=Man) 


Pulse  shaping  (l=y,  0=n) 


Num  taps  in  FIRs 


Num  taps  in  TIAF 


hJum  taps  in  each  bank  of  TDAF 


SOI  pulse  shaping  LPF  cutoff  freq 


SOI  pulse  shaping  LPF  gain 


SOI  carrier  amplitude 


SOI  carrier  frequency 


SNOI  baud  rale 


SNOI  carrier  amplitude 


SNOI  carrier  frequency 


outijut  LPF  gain 


output  LPF  cutoff  frequency 


demodulator  phase  shift 


Noi.se  gain 


misadjustment 


outputflag  (l=random  data.  0=square  \va.ve) 


Num  symbols  in  one  epoch 


Num  of  epochs  (does  not  include  adaptation) 


Number  of  adaptation  epochs 


Table  4.2.  Input  parameters  for  noisy  environment  simulations  for  the  BER.  version 


Input  Value 


6509731 


-1018 


64 


1 


1 


64 


32 


32 


1.0 


1.0 


See  Table 


6.0 


.95 


0.0 


6.0 


l.O 


1.0 
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Parameter 

6509731 

random  bit  generator  seed 

-1018 

AVVGN  seed 

64 

number  of  samples/symbol 

1 

Manchester  or  Bipolar  format  (l=Bi,  0=Man) 

1 

Pulse  shaping  (l=y,  0=n) 

64 

Num  taps  in  FIRs 

32 

Num  taps  in  TIAF 

32 

Num  taps  in  each  bank  of  TDAF 

1.0 

SOI  pulse  shaping  LPF  cutoff  freq 

1.0 

SOI  pulse  shaping  LPF  gain 

See  Table  4.I 

SOI  carrier  amplitude 

6.0 

SOI  carrier  frequency 

.95 

SNOI  baud  rate 

0.0 

SNOI  carrier  amplitude 

6.0 

SNOI  carrier  frequency 

1.0 

output  LPF  gain 

1.0 

output  LPF  cutoff  frequency 

0.0 

demodulator  phase  shift 

Sec  Table  JjA 

Noise  gain 

0.05 

misadjustment 

1 

oulputflag  {l=random  data,  0=squarc  wave) 

100 

Num  of  epochs 

16 

Number  of  symbols  to  average 

512 

Number  of  .symbols  per  epoch 
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Figure  4.1.  Improvement  factor  as  a  function  of  SNR 


4.2.1  Improvement  Factor.  Figure  4.1  shows  the  improvement  possible  for 
varying  SNR.  The  data  used  was  generated  by  the  BER  version  of  the  program  be- 
cuasc  that  version  freezes  the  filter  weights  after  a  specified  adaptation  time.  The 
data  for  Figure  4.1  was  taken  from  Table  4.4.  Reed  shows  in  his  PhD  dissertation 
(7:77)  that  the  TD.AF  (when  fulh'  adapted  -  what  he  calls  an  optimal  time  de¬ 
pendent  adaptive  filler)  will  always  do  at  le<ist  as  well  as  the  TIAF  (also  when  fully 
adapted).  Why,  then,  is  the  MSE  smaller  for  the  TI.4F  than  for  the  TDAF  at  20 
dB  in  Table  4.4?  The  key  is  that  tlie  TDAF  could  have  adapted  further  if  the  sim¬ 
ulation  had  not  frozen  its  filler  weights.  When  the  same  simulation  was  re-run  and 
the  adaptation  time  was  doubled  (from  512  adaptation  symbols  to  1024  adapation 
symbols),  the  MSE  for  the  TDAF  fell  to  0.002970,  and  the  improvement  factor  rose 
to  2.11212.  Therefore,  the  shape  of  the  curve  in  Figure  4.1  is  dependent  on  the 
amount  time  the  filter  is  allowed  to  adapt. 
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SNR 

MSE 

j 

TDAF. 

TIAF 

-20  dB 

0.550010 

-15  dB 

wimum 

0.484972 

mmii 

-10  dB 

0.272262 

Mmmm 

-5  dB 

0.141466 

■tiibtoaa 

0  dB 

0.064374 

0.125292 

1.94631 

5  dB 

Em 

0.060045 

2.20341 

10  dB 

0.012286 

0.026911 

2.19038 

15  dB 

0.007024 

0.011781 

1.67725 

20  dB 

0.005233 

0.86719 

Table  4.4.  MSE  data  for  simulations  in  AWGN 


The  symmetry  of  Figure  4.1  is  not  unexpected.  When  the  SNR  is  high,  both 
filters  achieve  a  nearly  optimum  solution,  resulting  in  an  improvement  factor  near 
unity.  When  the  SNR  is  low,  neither  filter  can  “lock  on”  to  the  SOI;  the  MSE 
for  both  filters  is  large,  and  again,  the  improvement  factor  drops  to  near  unity. 
When  the  SNR  is  in  the  moderate  range,  however,  the  TDAF’s  advantage  due  to 
stationarization  of  the  SOI  comes  into  play  as  illustrated  by  the  larger  improvement 
factor. 


4.2.S  Comparison  of  Filtered  Demodulated  Signals.  A  comparison  of  the  de¬ 
modulated  signals  for  various  levels  of  SNR  can  be  seen  in  Figures  4.2  through  4,6. 

Note  that  the  level  of  the  output  is  reduced  each  time  the  SNR  is  decreased.  The 
reason  for  that  is  simple.  In  order  to  keep  the  adaptation  time  the  same  for  each 
run,  the  total  amount  of  power  input  to  the  two  filters  is  normalized  to  0  dBW  (.sec 
Table  4.1).  Recalling  from  Section  2.3.2  that  the  adaptation  coefficient.  depends 
on  the  input  power  to  the  filter,  maintaining  an  input  power  of  0  dBW  to  the  filter 
results  in  a  constant  value  for  p  for  all  runs.  Hence,  the  adaptation  time  is  the  same. 

4.2.S  Comparison  of  MSE  Learning  Curves.  If  the  MSE  for  an  adaptive  filter 
is  plotted  as  a  function  of  sample  number,  tlie  result  is  the  “learning  curve”  for  the 
filter  (10:51).  For  the  TIAF,  the  data  for  the  learning  curve  is  collected  in  the 
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Volttt 


Figure  4.2.  The  desired  signal,  TDAF  response,  and  TIAF  response  for  0  dB  SNR 


Figure  4.-3.  The  desired  signal,  TDAF  response,  and  TIAF  response  for  -5  dB  SNR 
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Figure  4.6.  The  desired  signal.  TDAF  response,  and  TI.A.F  response  for  -20  dB  SNR 

following  manner.  The  filter  is  initialized  and  allowed  to  adapt:  the  resulting  error 
for  each  sample  is  squared  and  stored  in  an  array.  Then  the  filter  is  reinitialized,  and 
again  allowed  to  adapt.  The  resulting  error  is  squared,  and  added  to  the  previous 
error  data.  Aftei  this  process  is  repeated  M  limes,  where  M  is  some  sufficiently 
large  number  [M  =  100  in  this  thesis),  each  element  is  divided  by  M.  The  result  is 
an  e.stimate  of  the  MSE-.  The  larger  M.  the  better  the  estimate. 

For  Ihe  TDAF,  the  process  is  somewhat  more  complicated.  Recalling  from 
Section  2.3.3  and  3.3.7  that  the  individual  TIAFs  comprising  thcTD.AF  are  adapted 
only  once  per  symbol,  it  does  not  make  sense  to  plot  the  learning  curve  as  a  function 
of  absolute  sample  number.  This  is  because  no  filter  adaptation  has  occurred  from 
one  sample  to  the  next  in  the  TD.AF  as  is  docs  in  the  TIAF.  Instead,  the  adaptation 
for  the  TDAF  is  from  one  symbol  to  the  next.  Therefore,  for  the  TD.AF.  in  addition  to 
averaging  the  squared  error  for  each  reinitialized  run  to  the  next  as  in  th.e  TIAF,  the 
error  for  each  sample  comprising  a  given  symbol  is  averaged.  The  resulting  learning 
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Figure  4.7.  MSE  learning  curves  for  the  TDAF  and  TIAF  for  0  clB  SNR. 

curve  plots  the  MSE  as  a  function  of  symbol  number,  or  (sample  number)/ (samples 
per  symbol). 

Referring  to  Figures  4.7  through  4.11,  note  that  the  SNR  is  proportional  to 
the  initial  slope  of  the  learning  curve.  This  is  because  strong  correlation  between 
the  input  signal  and  the  desired  signal  yields  an  error  surface  with  steep  sides  (see 
Section  2.3.2).  Steep  error  surface  sides  result  in  faster  adaptation. 

The  learning  curve  plots  have  been  normalized  by  the  average  power  in  the 
desired  signal.  This  was  done  so  that  comparisons  of  the  rate  of  learning  could  be 
made  without  regard  to  the  magnitude  of  the  MSE. 

The  purpose  of  the  learning  curve  is  to  show  the  adaptation  process  of  the 
filters.  The  learning  curve  plots  are  not  intended  to  provide  an  indication  of  the 
value  of  the  MSE;  instead,  refer  to  Table  4.4. 


4-9 


Moan  SquAt'Atl  Ei'i'or 


Figure  4.8.  MSB  learning  curves  for  the  TDAF  and  TIAF  for  -5  clB  SNR 


Figure  4.9.  MSB  learning  curves  for  the  TDAF  and  TIAF  for  -10  dB  SNR 
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Figure  4.10.  MSB  learning  curves  for  the  TDAF  and  TIAF  for  -15  dB  SNR 


Figure  4.11.  MSB  learning  curves  for  the  TDAF  and  TIAF  for  -20  dB  SNR 


i 


SNR 

TDAF 

TIAF 

wmsm 

0.0 

0.0 

Bin 

0.004029 

0.008058 

-10  dB 

0.061653 

0.070321 

-15  dB 

0.198877 

-20  clB 

0.535764 

0.550010 

Table  4.5.  BER  for  five  simulations  in  a  noisy  environment 

J^.24  Bit  Error  Rate  An  error  detection  occurs  when  the  sum  of  samples 
over  a  symbol  time  is  less  than  zero  when  a  1  was  transmitted,  and  greater  than  zero 
when  a  -1  was  transmitted.  BER  is  calculated  according  to 


BER  = 


#  of  Bit  Errors 
#  of  Symbols  in  Simulation 


(4.2) 


Since  a  BPSK  modulatoi  was  used,  the  number  of  symbols  equals  the  number  of  bits. 
In  many  instances,  no  errors  were  detected  during  the  course  of  the  simulation.  It  is 
not  possible  to  overemphasize  the  point  that  the  BER  for  that  particular  SNR  is  not 
zero.  The  only  conclusion  that  may  be  made  is  that  for  that  particular  simulation 
run,  no  bit  errors  w’ere  made.  With  that  in  mind,  refer  to  Table  4.5.  Note  that 
the  greatest  reduction  in  BER  occurs  al  -5  dB  SNR.  These  results  seem  to  compare 
favorably  w'ith  those  plotted  in  Figure  4.1.  Each  simulation  w'as  8192  symbols  in 
duration. 


J,.2.5  Summarij  The  data  pre.sented  show-  that  the  TDAF  in  general  provides 
some  improvement  over  the  TI.AF  as  long  as  the  SNR.  is  no  worse  than  about  -10  dB. 
When  the  SNR  is  very  high,  both  filters  achie\e  a  low'  ,MSE;  hence,  the  improvement 
is  not  appreciable. 


4.S  Simulation  in  Interference 

In  this  section  the  performance  of  the  TDAF  is  measured  under  the  condition 
that  an  interference  signal  of  the  same  carrier  frequency  and  similar  baud  rate  is 
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SIR 

SOI 

-20  dB 

0.140720 

■gfiMiTtlil 

-15  dB 

0.247602 

1.392370 

-10  dB 

ligEaSTM 

1.348400 

-5  dB 

0.693186 

1.232678 

0  dB 

1.000000 

5  dB 

1.232678 

0.693186 

10  dB 

1.. 348400 

15  dB 

1.392370 

0.247602 

1.407195 

0.140720 

Table  4.6.  SOI  and  SNOI  carrier  amplitude  for  interference  environment 
simulations 

present  in  the  channel.  As  in  Section  4.2,  the  total  power  into  the  filter  is  maintained 
at  1  watt.  See  Table  4.6  for  the  carrier  amplitudes  used. 

Somewhat  unrealisticly,  the  condition  that  no  AWGN  is  present  on  the  channel 
was  chosen  in  an  attempt  to  identify  the  performance  differences  of  the  filters  for  noise 
and  interference.  The  ratio  of  interference  baud  rate  to  SOI  baud  rate  was  set  at  0.95 
for  all  simulations  in  this  section.  It  was  thought  that  such  a  ratio  would  adequately 
stress  the  filters  without  rendering  either  TDAF  or  TIAF  useless.  Surprisingly,  the 
adaptive  filters  are  relativeh'  insensitive  to  the  baud  rate  of  the  interfere!'  as  can  be 
seen  in  Figures  4.24  and  4.25.  Table  4.7  shows  the  input  parameters  used  to  obtain 
the  data  for  the  simulations  using  the  BER  version  of  the  program.  Table 4.8  shows 
the  input  parameters  used  for  the  LC  version  simulations. 

4-3.1  Improvement  Factor.  The  improvement  gain  for  the  TDAF  was  im¬ 
pressive  when  interference  was  jjreseni.  Table  4.9  shows  the  results  obtained  from 
the  BER  version  of  the  simulation.  The  data  shown  has  been  plotted  in  Figure  4.12. 

Note  the  similarity  in  shape  to  Figure  4.1  in  Section  4.2.1.  In  Figure  4.12  however, 
the  magnitude  of  the  improvement  is  significantly  higher  than  in  Figure  4.1.  This  is 
an  indication  that  the  major  advantage  of  the  TDAF  is  in  interference  rejection. 
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Input  Value 

Parameter 

6509731 

random  bit  generator  seed 

-1018 

AWGN  seed 

64 

number  of  samples/symbol 

1 

Manchester  or  Bipolar  format  (l=Bi,  0=Man) 

1 

Pulse  shaping  (l=y,  0=n) 

64 

Num  taps  in  FIRs 

32 

Num  taps  in  TIAF 

32 

Num  taps  in  each  bank  of  TDAF 

1.0 

SOI  pulse  shaping  LPF  cutoff  freq 

1.0 

SOI  pulse  shaping  LPF  gain 

See  Table  Jf.O 

SOI  carrier  amplitude 

6.0 

SOI  carrier  frequency 

.95 

SNOI  baud  rate 

See  Table  4.6 

SNOI  carrier  amplitude 

6.0 

SNOI  carrier  frequency 

1.0 

output  LPF  gain 

1.0 

output  LPF  cutoff  frequency 

0.0 

demodulator  phase  shift 

0.0 

Noise  gain 

0.05 

misadjustment 

1 

outputflag  (l=random  data,  0=square  wave) 

16 

Num  symbols  in  one  epoch 

512 

Num  of  epochs  (does  not  include  adaptation) 

32 

Number  of  adaptation  epochs 

b!e  4.7.  Input,  parameteis  lor  interference  environment  simulations  for  the  BER. 
version 
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Input  Value 

Parameter 

6509731 

random  bit  generator  seed 

-1018 

AWGN  seed 

64 

number  of  samples/symbol 

1 

Manchester  or  Bipolar  format  (l=Bi,  0=Man) 

1 

Pulse  shaping  (l=y,  0=n) 

64 

Num  taps  in  FIRs 

32 

Num  taps  in  TIAF 

32 

Num  taps  in  each  bank  of  TDAF 

1.0 

SOI  pulse  shaping  LPF  cutoff  freq 

1.0 

SOI  pulse  shaping  LPF  gain 

See  Table  4.6 

SOI  carrier  amplitude 

6.0 

SOI  carrier  frequency 

.95 

SNOI  baud  rate 

See  Table  4.6 

SNOI  carrier  amplitude 

6.0 

SNOI  carrier  frequency 

1.0 

output  LPF  gain 

1.0 

output  LPF  cutoff  frequency 

0.0 

demodulator  phase  shift 

0.0 

Noise  gain 

0.05 

misadjustment 

1 

outputflag  (l=random  data,  0=square  wave) 

100 

Num  of  epochs 

16 

Number  of  symbols  to  average 

512 

Number  of  symbols  per  epoch 

Table  4.S.  Input  parameters  for  interference  environment  simulations  for  the  LC 
version 


SNR 

MSE 

.1 

TDAF 

TIAF 

-20  dB 

0.455682 

0.954006 

2.09358 

-15  dB 

0.269755 

0.870597 

3.22736 

-10  dB 

0.096536 

0.741551 

7.68160 

-5  dB 

0.0.39515 

0.570294 

14.4.3234 

OdB 

0.019807 

0.370759 

18.71878 

5dB 

0.011584 

0.183849 

15.87094 

10  dB 

0.008921 

0.072100 

8.08285 

15  dB 

0.009487 

0.028873 

3.04343 

20  dB 

0.007489 

0.012307 

1.64334 

Table  4.9.  MSE  data  for  simulations  in  interference 
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Figure  4.12.  Improvement  factor  as  a  function  of  SIR 


4.3.2  Comparison  of  filtered  Demodulated  Signals.  A  comparison  of  the  cle- 
moclulatecl  signals  for  various  levels  of  SIR  can  be  seen  in  Figures  4.13  through  4.17. 
Note  that  as  the  SIR  gets  smaller,  the  response  of  the  TIAF  loses  all  similarity  to  the 
desired  signal.  Foi  the  TDAF,  on  the  other  hand,  the  shape  of  the  response  is  very 
nearly  correct  even  at  -20  dB  SIR,  although  the  level  of  its  output  is  not.  Figure  4.18 
is  a  plot  of  the  SNOI  over  the  same  time  interval  as  Figures  4.13  through  4.17.  Note 
the  similarity  of  the  TI.AF  output  in  Figures  4.16  and  4.17  to  the  plot  in  Figure  4.18 
indicating  that  the  TI.AF  could  not  track  tlie  SOI,  and  instead  pa,ssed  the  SNOI. 

4.3.3  Comparison  oj  MSB  Learning  Curves.  The  learning  curves  shown  in 
Figures  4.19  through  4.23  arc  computed  in  the  same  fashion  as  those  in  Section  4.2.3. 
Again,  the  adaptive  filters  converge  more  rapidly  when  the  SIR  is  high.  In  fact 
Figures  4.22  and  4.23  show  that  there  is  little  reduction  in  the  level  of  the  MSE 
throughout  the  course  of  the  simulation. 
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Figure  4.13.  The  desired  signal.  TDAF  response,  and  TIAF  response  for  0  dB  SIR 


Figure  4.14.  The  desired  signal,  TD.AF  response,  and  TIAF  response  for  -5  dB  SIR 


Figure  4.15.  The  desired  signal,  TDAF  response,  and  TIAF  response  for  -10  dB 
SIR 
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Figure  4.17.  The  desired  signal,  TDAF  response,  and  TIAF  response  for  -20  dB 

SIR 


Figure  4. IS.  The  SNOI  present  in  Figures  4.13  through  4.17 


Figure  4.20.  MSB  learning  curves  for  the  TD.^VF  and  TIAF  for  -5  dB  SIR 
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Figure  4.21.  MSE  learning  curves  for  the  TDAF  and  TIAF  for  -10  dB  SIR 


Figure  4.22.  MSE  learning  curves  for  the  TD.AF  and  TI.AF  for  -15  dB  SIR 
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SIR 

TD.AF 

TIAF 

OdB 

0.0 

0.220974 

-5  dB 

0.0 

0.396044 

-10  dB 

0.0 

0.449152 

-15  dB 

0.0 

0.46770S 

-20  clB 

0.152240 

0.481016 

Tabic  4.10.  BER  for  five  Simulations  in  Interference 


4-3.4  Bil  Error  Rate  Table  4.10  shows  the  BER  observed  during  the  simula¬ 
tions  in  interference.  These  five  simulations  indicate  that  there  may  be  a  significant 
improvement  in  BER  when  a  TD.AF  is  used  to  e.vtract  a  signal  from  a  strong  inter- 
ferer  .  Referring  to  Table  I..55  note  that  in  noise,  both  filters  suffered  approximately 
50%  error  rate.  In  interference,  however,  the  BER  for  the  TDAF  dropped  from 
appro.xi match'  50%  to  under  20%. 


4.3.0  Summary.  The  improvement  obtained  for  simulations  in  an  interfer¬ 
ence  were  better  overall  than  for  those  run  in  a  noisy  environment.  The  improvement 
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factor  was  higher  at  all  levels  of  SIR,  and  the  Bit  Error  Rate  appeared  to  be  lower. 
The  simulations  show  that  the  TDAF  is  very  good  at  rejecting  interference,  even 
when  the  interfere!-  is  at  the  same  carrier  frequency. 

jfJf  Varyi.'j  the  Baud  Rate  of  the  Inierferer 

The  data  used  to  produce  Figure  4.24  and  4.25  was  compiled  from  12  simula¬ 
tions.  A  simulation  was  run  for  each  the  four  relative  baud  rates  at  -10  dB,  0  dB  and 
10  dB  SIR.  The  carrier  amplitudes  for  those  levels  of  SIR  can  be  found  in  Table  4.6. 
All  remaining  input  parameters  were  the  same  as  those  listed  in  Table  4.7. 

The  data  in  this  section  support  the  surprising  result  that  the  TDAF  is  es¬ 
sentially  unaffected  by  the  baud  rate  of  the  interfere!-.  Figure  4.24  illustrates  the 
TDAF’s  ability  to  reject  an  SNOI  even  when  its  baud  rate,  modulation  type  and 
carrier  frequency  are  the  same  as  the  SOI.  Figure  4.25  is  included  to  show  the  TIAF’s 
performance  under  the  same  conditions.  The  TIAF  was  also  fairly  insensitive  to  the 
baud  of  the  interfere!-,  but  note  that  the  MSE  for  the  TIAF  is  approximately  10  dB 
worse. 


f.S  Varying  the  Carrier  Frequency  of  the  SNOI 

The  data  for  Figure  4.26  were  compiled  from  five  simulations.  In  each  simula¬ 
tion,  the  carrier  frequency  of  the  SNOI  was  varied  as  indicated.  The  SIR  was  fixed 
at  0  dB.  The  '-emaining  input  parameters  to  the  simulation  were  set  according  to 
Table  4.7. 

The  TDAF  is  not  affected  by  SNOI  carrier  frequency  proximity  as  illustrated 
by  Figure  4.26.  The  TIAF  however  suffers  a  significant  increase  in  MSE  as  the  carrier- 
frequencies  of  the  SOI  and  SNOI  coincide. 
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Figure  4.24.  MSE  as  a  function  relative  baud  rate  at  various  SIRs  for  the  TDAF 


Figure  4.25.  MSE  as  a  function  relative  baud  rate  at  various  SIRs  for  the  TIAF 


Figure  4.26.  MSE  as  a  function  SNOI  carrier  frequency.  SOI  carrier  frequency  fixed 
at  6  Hz 

4.6  Chapter  Summary 

Over  50  simulations  were  run  to  obtain  the  data  for  this  chapter.  The  results 
of  those  simulations  indicate  that  that  for  a  given  communicatiojis  system,  if  a  choice 
exists  between  a  TDAF  and  a  TIAF,  then  the  conditions  under  which  the  filter  will 
be  operating  must  be  considered  prior  to  making  the  selection.  The  simulations 
yielded  improvement  factors  that  ranged  from  -0.6  dB  to  12  dB.  Even  though  there 
was  one  case  where  the  TIAF  outperformed  the  TDAF,  it  was  shown  that  given 
sufficient  adaptation  time,  the  TD.'\F  will  always  do  at  least  as  well  as  the  T1.4F. 


4-25 


V.  Conclusions  and  Recommendations 


5. 1  Conclusions 

This  thesis  presented  a  TDAF  which  could  be  used  to  improve  the  signal 
to  interference  ratio  (SIR)  and  signal  to  noise  ratio  (SNR)  of  digitally  modulated 
communications  signals.  The  performance  improvement  of  the  TDAF  over  the  TIA? 
was  determined  based  on  the  application  of  various  metrics,  including  Mean  Square 
Error  (MSE)  and  Bit  Error  Rate  (BER).  The  collected  data  indicate  that  there  are 
advantages  and  disadvantages  of  the  TDAF.  The  advantages  can  be  summarized  as 
follows: 

•  The  improvement  in  MSE  obtained  by  using  a  TDAF  over  a  TIAF  in  a  noisy 
environment  can  be  significant.  With  a  0  dB  SNR,  a  performance  gain  of 
nearly  3  dB  can  be  expected. 

•  The  improvement  in  MSE  obtained  by  using  a  TD.AF  over  a  TIAF  in  situation 
where  there  is  a  strong  interference  signal  can  be  even  more  significant.  With 
a  0  dB  SIR,  a  performance  gain  of  more  than  12  dB  can  be  expected. 

•  A  receiver  fitted  with  a  TD.A.F  which  has  been  given  adequate  time  to  adapt 
is  likely  to  have  a  lower  BER  than  a  receiver  using  a  TLAF. 

•  The  TDAF  always  does  at  least  as  well  as  the  TIAF,  given  adequate  time  to 
adapt. 


The  disadvantages  are: 

•  The  TDAF  takes  consicloiably  longer  to  adapt  than  the  TIAF.  While  the  TIAF 
adapts  at  each  sample  time,  the  TDAF  actually  only  “fractionally"  adapts  at 
each  sample  time.  A  full  symbol  is  required  to  adapt  each  TIAF  in  the  TD.AF 
structure. 
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•  The  TDAF  would  require  much  more  “real  estate”  if  implemented  in  silicon. 
This  is  due  to  the  multiplicitj'  of  filter  coefficients  required  for  the  TDAF. 
Fortunately,  there  is  little  added  computational  overhead,  since  the  process  of 
updating  filter  coefficients  is  no  more  frequent  for  the  TDAF  than  it  is  for  the 
TIAF. 

•  The  TDAF  is  somewhat  more  complex  than  the  TIAF.  This  should  be  clear 
since  the  TDAF  is  made  up  of  a.  bank  of  parallel  TIAFs. 

•  As  SNR  at  the  filtei  input  decreases,  the  improvement  in  Mean  Square  Error 
decreases.  By  the  time  the  SNR.  has  degraded  to  less  than  -10  dB,  the  per¬ 
formance  gain  has  dropped  to  under  1  dB.  If  that  is  the  case,  then  given  the 
other  disadvantages  of  the  TDAF,  the  use  of  a  TIAF  might  make  more  sense. 

5.%  Recommendations 

There  are  several  topics  that  could  (and  probably  should)  have  been  covered 
in  the  execution  of  this  research.  Some  of  the  more  urgent  are  listed  below. 

1.  Investigate  the  role  that  the  modulation  type  plays  in  the  performance  gain. 
By  implementing  different  modulator/demodulator  pairs,  the  simulation  can 
be  used  to  determine  im])rovement  gain  as  a  function  of  modulation  type. 
Other  modulation  schemes  might  include: 

•  Quadrature  Pha.se  Shift  Keying  (QPSK) 

•  Offset  Quadrature  Phase  Shift  Keying  (OQPSK) 

•  M-ary  Frequency  Shift  Keying  (FSK) 

•  Minimum  Shift  Keying  (MSK) 

•  Quadrature  Amplitude  Modulation  (QAM) 
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2.  Implement  the  TDAF  with  no  external  training  signal.  With  some  added 
complexity,  the  inconvenience  of  having  to  transmit  a  known  preample  can  be 
circumvented  (7),  and  the  filter  can  be  allowed  to  adapt  continuously. 

3.  Implement  the  TDAF  in  the  frequency  domain.  The  advantage  of  doing  so  can 
be  significant.  The  cycle  frequencies  of  the  signal  used  in  this  simulation  were 
at  ±2 fc-  In  the  frequency  domain,  the  optimum  TDAF  is  implemented  with  as 
many  TIAFs  as  exist  cycle  frequencies  (7:103).  One  advantage  of  fewer  TIAFs 
may  be  more  rapid  adaptation  since  there  would  be  fewer  weight  vectors  to 
update. 

4.  Precalculate  an  initial  weight  vector  for  the  adaptive  filters.  Both  TDAF  and 
TI.AF  are  initialized  with  zero  weight  vectors  in  the  simulation  presented  here. 
It  should  be  possible  to  calculate  a  “near  optimum”  solution  for  the  adaptive 
filters  for  a  given  SOI,  and  then  allow  the  filter  weights  adapt  from  that  point. 

5.  Investigate  the  apparent  inconsistency  that  exists  between  the  MSB  returned 
for  the  LC  version  and  the  MSB  returned  for  the  BBR  version.  While  the  BER 
consistantly  returns  smaller  Bit  Error  Rate  numbers  and  MSB  for  the  TDAF, 
the  -MSB  returned  by  the  LC  version  is  considerably  higher  for  the  TDAF. 

6.  Determine  the  improvement  of  the  TDAF  over  the  TIAF  for  the  important 
case  of  multiple  interferers. 

7.  Determine  the  improvement  of  the  TDAF  over  the  TI.AF  under  the  conditions 
of  simultaneous  interference  and  noise. 
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Appendix  A.  Inpvi  Parameters  and  Output  Files  for  the  BER 

Version 

The  BER  version  of  the  simulation  prompts  the  user  for  the  following  input 
parameters  at  run  time; 

•  Randovi  bit  generator  seed.  Seeds  the  random  bit  generator. 

•  Random  number  generator  seed.  Seeds  the  AVVGN  generator. 

•  Data  Sample  Rate.  The  number  of  samples  per  bit  of  the  data. 

•  Format.  \  flag  to  select  between  bi-polar  and  bi-phase  baseband  format. 

•  Pulse  shaping.  A  flag  to  allow  or  disallow  bandlimiting  of  the  baseband  data 
signal. 

•  FIR  taps  The  number  of  coefficients  in  tiie  non-adaptive  FIR  filters.  All  FIR 
filters  in  the  simulation  then  have  the  same  number  of  taps. 

•  TIAF  taps.  The  number  of  filler  coefficients  for  the  time  indejjendent  adaptive 
filter.  Stated  another  way,  the  number  of  elements  of  W. 

»  TD/\F  taps.  The  number  of  filter  coefficients  for  the  time  dependent  adaptive 
filter.  See  Section  2.3.3 

•  Bandlimiting  LPF  cutoff  frequency.  Passed  to  the  simulation  in  hertz.  Both 
SOI  and  SNOI  channels  use  the  same  cutoff. 

•  Bandlimiting  LPF  gain.  Allows  adjustment  of  the  level  of  the  output  of  the 
pulse  shapers. 

•  SOI  carrier  amplitude. 

•  SOI  carrier  frequency. 
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•  SNOI  symbol  frequency.  Expressed  as  a  fraction  of  the  SOI  symbol  frequency 
of  1  symbol  per  second.  May  be  any  positive  real  number. 

•  SNOI  carrier  amplUude. 

•  SNOI  carrier  frequency. 

•  Output  LPF  gain.  Allows  adjustment  of  the  demodulator  output  level. 

•  Output  LPF  cutoff  frequency.  In  hertz. 

•  Demodulator  phase  shift.  .Allows  for  adjustment  of  the  phase  of  the  bandpass 
signal  into  the  demodulator. 

•  Noise  gain.  Allows  for  adjustment  of  noise  power  in  the  unfiltered  bandpass 
signal. 

•  Misadjustment.  See  Section  2.3.3 

•  Data  type  fl.  q.  Allows  for  selection  between  a  square  wave  baseband  signal  or 
random  data. 

•  Symbols  per  epoch.  The  program  writes  this  number  of  .symbols  to  a  data  file 
at  tlie  end  of  a  simulation  run. 

•  Ntimber  of  run  epochs.  Used  to  specify  the  number  of  .symbols  for  the  entire 
simulation  (total  symbols  =  symbols  per  epoch  x  number  of  run  ep'^chs). 

•  Number  of  adaptation  epochs.  Used  to  specify  the  length  of  time  the  filters  arc 
allowed  to  adapt. 

Each  run  of  the  BER.  Version  produces  six  data  files.  The  files  are  written  into 
the  current  directory.  They  are: 

•  tdafvoll.dat  The  output  of  the  TDAF  over  the  last  specified  number  of  .symbols 
(see  above). 
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tiafvolLdat  The  output  of  the  TIAF  over  the  same  period. 

dsrdvoll.dat  The  desired  filter  output  over  the  same  period. 

tdafvec.dat  The  error  signal  of  the  TDAF  over  the  same  period. 

liafvic.dat  The  error  signal  of  the  TIAF  over  the  same  period. 

nunibers.lex  A  file  that  contains  various  data  and  figures  of  merit  for  the 
simulation  run: 

1.  The  adaptation  coefficient  for  the  simulation  (see  Section  2. -3. 3). 

2.  The  number  of  adaptation  epochs  in  the  simulation. 

3.  The  total  number  of  epochs  in  the  simulation. 

4.  The  number  of  post-adaptation  symbols  in  the  simulation. 

5.  The  mean  squared  error  for  the  TDAF. 

6.  The  mean  squared  error  for  the  TIAF. 

7.  The  number  of  symbol  errors  for  the  TDAF. 

S.  The  number  of  symbol  errors  for  the  TIAF. 


Appendix  B.  Input  Parameters  and  Output  Files  for  the  LC  Version 


The  LC  version  of  the  simulation  prompts  the  user  for  the  following  input 
parameters  at  run  time: 

•  Random  bit  generator  seed.  Seeds  the  random  bit  generator. 

•  Random  number  generator  seed.  Seeds  the  AVVGN  generator, 

•  Data  .Sample  Rate.  The  number  of  samples  per  bit  of  the  data. 

•  Formal.  A  flag  to  select  between  bi-polar  and  bi-phase  baseband  format. 

•  Pulse  shaping.  A  flag  to  allow  or  disallow  bandlimiting  of  the  beiseband  data 
signal. 

•  FIR  taps  The  number  of  coefficients  in  the  non-adaptive  FIR  filters.  .411  FIR 
filters  in  the  simulation  then  have  the  same  number  of  taps. 

•  TI.AF  taps.  The  number  of  filter  coefficients  for  the  time  independent  adaptive 
filter.  Stated  another  way,  the  number  of  elements  of  W. 

•  TDAF  laps.  The  number  of  filter  coefficients  for  the  time  dependent  adaptive 
filter.  See  Section  2.3.3 

•  Bandlimiting  LPF  cutoff  frequeneg.  Passed  to  the  simulation  in  hertz.  Both 
SOI  and  SNOl  channels  use  the  same  cutoff. 

•  Bandlimiting  LPF  gain.  .Allows  adjustment  of  the  level  of  the  output  of  the 
pulse  shapers. 

•  .SOI  carrier  amplitude. 


•  SOI  can'ier  frequency. 


•  SNOl  symbol  frequency.  Expressed  as  a  fraction  of  the  SOI  syrnbol  frequency 
of  1  symbol  per  second.  May  be  any  positive  real  number. 

•  SNOI  caiTier  mnplitude. 

•  SNOI  cairier  frequency. 

•  Output  LPF  gain.  Allows  adjustment  of  the  demodulator  output  level. 

•  Output  LPF  cutoff  frequency.  In  hertz. 

•  Demodulator  phase  shift.  Allows  for  adjustment  of  the  phase  of  the  bandpass 
signal  into  the  demodulator. 

•  Noise  gain.  .Allows  for  adjustment  of  noise  power  in  the  unfiltered  bandpass 
signal. 

•  Misadjustmenl.  See  Section  2.3.3 

•  Data  type  flag.  Allows  for  selection  between  a  square  wave  baseband  signal  or 
random  data. 

•  Number  of  epochs.  The  greater  the  number  of  epochs,  the  better  the  estimate 
of  the  MSE. 

•  Number  of  symbols  to  average.  The  MSE  from  this  number  of  symbols  will  be 
averaged  at  the  end  of  each  epoch  to  provide  an  estimate  of  the  MSE. 

•  .Number  symbols  per  epoch.  This  number  of  symbols  will  be  applied  to  the 
TDAF  and  TI.AF  each  epoch. 

Each  run  of  version  LC  produces  three  data  files.  The  files  are  written  into  the 
current  director}'.  They  are: 

•  tdaflm.dat  An  estimate  of  the  expected  value  of  the  MSE  as  a  function  of  time 
for  the  TDAF. 
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tiajlrn.dat  An  estimate  of  the  expected  value  of  the  MSE  as  a  function  of4ime 
for  the  TIAF. 

hivmbers.lex  A  file  that  contains  various  data  and  figures  of  merit  for  the 
simulation  run: 

1.  The  adaptation  coefficient. 

2.  The  number  of  sj'mbols  represented  by  the  adaptation  period. 

3.  The  mean  squared  error  of  the  TDAF  for  the  last  n  symbols  of  the  epoch 
(where  n  is  selectable  at  run  time). 

4.  The  mean  squared  error  of  the  TIAF  for  the  last  n  symbols  of  the  epoch 
(where  n  is  selectable  at  run  time). 


Appendix  C.  Source  Code  for  the  BER  Version 


lie  :(c:jc;|c*:|c*;|c******«**^******4<************«***«*=l‘*3|‘***«*««****«**4:**/ 

/******  This  program  is  a  simulated  digital  communication  system  ******/ 

/* 

/  5  Oct  91:  This  version  of  the  program  has  an  optimized  TDAF. 

/  The  filter  now  has  a  one  dimensional  vector  (rather  than  a  matrix 

/  as  in  earlier  versions)  to  hold  the  values  of  the  input.  I  made 

/  the  startling  (to  me)  observation  that  each  row  of  the  input 
/  matrix  was  identical. 

/  7  Oct  91:  Added  input  parameter  "frozen"  to  the  adaptive 

/  filters.  When  frozen  =  0,  weights  are  allowed  to  adapt. 

/  After  the  filters  have  been  given  the  specified  length  of  time 
/  to  adapt,  then  frozen  is  set  to  1,  and  filter  weights  are  no 
/  longer  allowed  to  adapt. 

/  Executable:  frzloop 

/  */ 
/#*****+*♦*!(:  *****!►***************♦:***♦*****♦♦*♦*♦**********!(=******!(=♦♦**/ 


#define  spare  1 

^include  <stdio.h> 

S include  <math.h> 
^include  <stdlib.h> 

^define  pi 
#define  data_freq 
#define  data_rate 
Sdefine  sqr(x) 


/*  *f 


(4*atan(l .0) ) 
1.0 

1.0/datafreq 

(x)*(x) 


/:(:)(:*****:(!*:*::*:*♦*******♦**  FunCtiOU  PrOtOtypeS  *1F^^*****W*i^^!^***'*****^^^•****■/ 

double  setcutoff (double  samplerate, 
double  cutoff) ; 

void  calcfilterweights(int  numdelays, 

double  oraega.c. 


double  ^contents, 
double  ^weights) ; 

int  irbit2 (unsigned  long  *bseed) ; 

int  datagen(double  time, 

double  datarate, 

int  lastdata, 

double  *lasttime, 

unsigned  long  *bseed, 
int  outputf lag) ; 

double  Ipf (double  input, 

int  numdelays , 

double  ^weights , 

double  *contents , 

double  gain) ; 

double  modulate (double  input, 

double  carrierfreq, 
double  carrierampl , 
double  time) ; 

void  inittiaf (double  *weights, 
double  *contents, 
int  numtaps , 
double  terror, 
double  *bias) ; 

double  tiaf (float  input, 

float  desired, 

float  mu, 

int  numtaps , 

double  *error, 

double  *weights, 

double  ^contents,  double  *bias, 

int  frozen) ; 

double  demodulate (double  input. 

double  carrierfreq, 
double  phase, 
double  time) ; 


float  ranl(int  *idura); 


float  gasdev(int  *idum) ; 

double  **inittdafwts(int  sampersym, 
int  numtaps) ; 

double  *inittdaf contsdnt  numtaps); 

double  tdaf (double  input, 

double  desired, 
double  mu, 
int  sampersym, 
int  numtaps , 
double  **tdafwts, 
double  *tdaf cents, 
double  terror, 
double  *bias, 
int  frozen) ; 

double  *initerror(int  sampersys); 

double  *initbias(int  sampersym); 

void  nrerror(char  error.text  [] ) ; 

double  *dvector(int  nl, 
int  nh) ; 

double  **dmatrix(int  nrl, 
int  nrh, 
int  ncl, 
int  nch) ; 

float  *vector(int  nl, 
int  nh) ; 

int  manchesterCint  input, 

double  time, 
double  datarate, 
double  lasttime); 

void  free„dvector(double  *v, 
int  nl, 
int  nh) ; 

void  free_dmatrix (double  **m, 
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int  nrl, 
int  nrh, 
int  ncl, 
int  nch) ; 

int  bipolarCint  input); 


int  mainO  { 


unsigned  long 
numsamples , 
bbeed, 
bitvar=0; 

int  numsymbols, 
numf irtaps , 
numtiaftaps, 
numtdaftaps, 
numf irdelays , 
numtiaf delays , 
numloops , 
samplefreq, 
count , 
outputflag, 
adaploops , 
sample, 
s ample 1 , 
loopnum, 
dlt=l , 
d2t, 
d3t=l, 
d4t , 
d5t , 
d6t , 

sampersym, 
format , 
shape, 

tdaf biterror=0 , 
tiafbiterror=0, 
offsetfound=0, 
xsample. 
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bitcount  =  0, 
avgloops , 
msecount  =  0, 
nseed , 
'frozen=0: 


double 

samplerate , 

datafreq  =  data.freq, 

datar'ate  =  data.rate, 

carrierfreq, 

carrierarapl , 

noisegain, 

misadjust , 

mu , 

deltaphase , 
gainl, 

_  gain2, 
tiafgain2, 
avg  =  0.0, 
tiafavg  =  0.0, 

Ipf in.cutoff , 

Ipf out.cutof f , 
ilpfin.cutoff , 

*lpf conts , 

♦Ipfwts , 

♦tiafwts , 

’•‘tiaf  conts , 

♦nnwts , 

♦nnconts , 

♦  outweights , 
★outcontents , 

*tiaf outweights , 

♦tiaf outcontents , 
tiaf error  =  0.0, 
fiafbias  =  0.0, 
omega.c , 
omega.cl , 
time  =  0.0, 
sit,  s2t,  s3t, 
s4t,s5t,  s6t, 
s7t,  s8t , 
x6t,  x7t,  x8t, 
nnlt,  nn2t, 
tdafbitnum=0.0, 
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tiafbitnuin=0.0, 
dsrdbitnuin=0.0, 
s3tlast , 
nn2tlast, 
lasttime, 

**tdafwts, 

*tdaf cents , 

*tdaferror , 

*tdafbias, 

idatarate, 

idatafreq, 

ilasttime, 

*ilpf cents,  *ilpfwts, 
igainl , 
icarrierfreq, 
icarrierampl, 
itime  =  0.25, 
iemega_c, 
desiredpewer, 
interfpewer ; 

fleat 

*mse,  *tempmse, 

*t i af ms  e ,  *t i af t  empms  e , 
adapf acter=0 . 1 178 , 
leopfacter ; 

char 

buffer[128] ; 

FILE 

*tdafvelt,  *tiafvelt, 

*dsrdvclt,  *tdafvec, 

*tiafvec,  ^numbers; 

if  ((tdafvelt  =  fepenC'tdaf velt .dat" ,  "w"))  ==  NULL) 
printf  ("  ***  Ceuld  net  epen  tdafvelt.dat!  ***  \n"); 

if  ((tiafvelt  =  fopenC'tiaf velt. dat",  "w"))  ==  NULL) 
printf  ("  ***  Ceuld  net  epen  tiafvelt.dat!  ***  \n"): 

if  ((dsrdvelt  =  fepenC'dsrdvelt.dat",  "w"))  ==  NULL) 
printf  ("  ***  Could  not  open  dsrd/olt.dat!  ***  \n"); 

if  ((tdafvec  =  fopen("tdafvec.dat" ,  "w"))  ==  NULL) 
printf  ("  ***  Could  not  open  tdafvec.dat!  ***  \n") ; 
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if  ((tiafvec  =  fopenC'tiafvec.dat" ,  "w"))  ==  NULL) 
printf  ("  ***  Could  not  open  tiafvec.dat!  ***  \n"); 

if  ((numbers  =  f open ( "numbers. tex" ,  "w"))  ==  NULL) 
printf  ("  ***  Could  not  open  tdafnum.dat!  ***  \n"); 


User  Input  Section 

printf ("Seed  the  random  bit  generator:  "); 
gets (buffer) ; 

sscanf  (buffer,  "*/,U",  &bseed)  ; 
printf  ("*/,d\n" ,  bseed) ; 

printf ("Seed  the  AWGN  Generator  (integer  <  0):  "); 
gets (buffer) ; 

sscanf  (buffer,  "'/,d",  fenseed)  ; 
printf  ("y,d\n" ,  nseed) ; 

printf ("Number  of  samples  per  symbol  required:  "); 
gets (buffer) ; 

sscanf  (buffer,  "'/,d",  fesampersym); 
printf  ("7, d\n" ,  sampersym) ; 

samplefreq  =  sampersym*datafreq; 
samplerate  =  (double)  1 .0/samplefreq; 

printf ..unchester  (  ^  ■>.  Bipolar  (1)  format:  "); 
gets (buffer) ; 

sscanf  (buffer,  "7,d",  ftformat); 
printf  ("V.dXn" ,  format) ; 

printf ("Pulse  shaping?  (l=y,  0=n):  "); 
gets (buffer) ; 

sscanf  (buffer,  "7,d",  &shape) ; 
printf  ("V.dXn" ,  shape); 

printf ("Number  of  taps  in  the  FIR  filters:  "); 
gets (buffer) ; 

sscanf  (buffer ,  "7.d",  &numf  irtaps) ; 
printf  ("V.dXn" ,  numf  irtaps) ; 


numfirdelays  =  numf irtaps  -  1; 


printf ("Number  of  taps  in  the  TIAF  adaptive  filter:  "); 
gets (buffer) ; 

sscanf  (buffer,  "*/,d",  ftnumtiaftaps) ; 
printf  ("*/,d\n" ,  numtiaftaps) ; 

numtiafdelays  =  numtiaftaps-1 ; 

printf ("Number  of  taps  in  the  TDAF  adaptive  filter:  "); 
gets (buffer) ; 

sscanf  (buffer,  "*/,d",  &numtdaftaps) ; 
printf  ("'/,d\n" ,  numtdaftaps) ; 

printf ("SOI  pulse  shaping  LPF  cutoff  (Hz):  "); 
gets (buffer) ; 

sscanf  (buffer ,  '“/.If",  &lpfin_cutoff  )  ; 
printf  ("*/,f\n" ,  Ipf  in_cutoff ) ; 

ilpfin_cutoff  =  Ipfin.cutoff ; 

printf ("Gain  of  pulse  shaping  LPF:  "); 
gets(buff er) ; 

sscanf  (buffer,  "’/.If",  &gainl); 
printf  ("7.f\n",  gainl); 

igainl  =  gainl; 

printf ("SOI  carrier  amplitude:  "); 
gets (buffer) ; 

sscanf  (buffer ,  "'/.If",  ftcarrierampl); 
printf  ("'/,f\n"  ,  carrierampl) ; 

desiredpower  =  (sqr(carrierampl))/2; 

printf ("SOI  carrier  frequency:  "); 
gets(buff er) ; 

sscanf  (buffer ,  "'/.If",  fccarrierfreq) ; 
printf  ("'/,f\n" ,  carrierfreq) ; 

printf ("Symbol  frequency  for  the  interferer:  "); 
gets (buffer) ; 

sscanf  (buffer,  "’/.If",  feidatafreq) ; 
printf  ("'/,f\n"  ,idatafr3q) ; 

idatarate  =  1 .0/idatafreq; 
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printfC'SNOI  carrier  amplitude:  "); 
gets (buffer) ; 

sscanf  (buff  er,  '"/.If",  fticarrierampl) ; 
printf  ("7,f\n" ,  icarrierampl) ; 

interf power  =  (sqr( icarrierampl) )/2; 

printfC'SNOI  carrier  frequency:  "); 
gets (buffer) ; 

sscanf (buffer,  "7lf",  fticarrierfreq) ; 
printf  ("7, f\n" ,  icarrierfreq)  ; 

printf ("Gain  of  output  LPF:  ") ; 
gets (buffer) ; 

sscanf  (buffer,  "7.1f",  &gain2) ; 
printf  ("7.f\n" ,  gain2) ; 

printf ("Output  LPF  cutoff  (Hz):  "); 
gets (buffer) ; 

sscanf  (buffer,  "7ilf",  &lpf  out_cutoff  )  ; 
printf  ("7if\n" ,  lpfout_cutoff ) ; 

printf ("Phase  shift  for  demodulator:  "); 
gets (buffer) ; 

sscanf  (buffer ,  "7.1f",  &deltaphase)  ; 

printf  ("7.f\n",  deltaphase);  , 

printf ("Noise  factor:  "); 
gets(buff er) ; 

sscanf  (buffer,  "7.1f",  &noisegain)  ; 
printf  ("7.f\n",  noisegain)  ; 

printf ("Misadjustment  factor:  "); 
gets (buffer) ; 

sscanf  (buffer,  "7ilf",  ftmisadjust) ; 
printf ("7.f\n" ,  misadjust); 

if  (carrierfreq  >  0.0) 

mu  =  misadjust/((desiredpower  +  interf power  + 
sqr (noisegain) )*(numtdaf taps)) ; 

else 

mu  =  misadjust/ ((sqr(carrierampl)  +  sqr (icarrierampl) 
sqr(noisegain))*(numtdaftaps)) ; 


fprintf (numbers , "Mu  =  '/,f\n",  mu); 

printf ("Random  data  (1)  or  square  wave  (0):  "); 

gets(buff er) ; 

sscanf  (buffer,  "'/,d",  ftoutputf lag)  ; 
printf  ("7.d\n" ,  outputf  lag) ; 

printf ("Number  of  symbols  in  one  epoch:  "); 
gets (buffer) ; 

sscanf  (buffer,  '"/.u",  fenumsymbols)  ; 
printf  ("'/,d\n" ,  numsymbols) ; 

numsamples  =  numsymbols+sampersym; 

adaploops  =  ((adapf actor /mu) /numsymbols) ; 

printf ("Number  of  epochs  to  average:  "); 
gets (buffer) ; 

sscanf  (buffer,  '"/.d"  ,  feavgloops); 
printf  ("‘/.d\n" ,  avgloops)  ; 

printf ("Number  of  adapation  epochs:  ") ; 
gets (buffer) ; 

sscanf  (buffer,  "7,d",  feadaploops)  ; 
printf  ("7.d\n" ,  adaploops) ; 

numloops  =  adaploops  +  avgloops; 

*printf (numbers , "Number  of  adaptation  epochs:  7.d\n",  adaploops); 
fprintf  (numbers, "Total  number  of  epochs:  7.d\n",  numloops); 
printf  ("Total  number  of  epochs:  7,d\n",  numloops); 
fprintf  (numbers,  "Total  number  of  post-adaptation  symbols:  7id\n", 
numsymbols*avgloops) ; 
loopf actor  =  1 .©/(numloops) ; 

Initialization  Section 

************************■********■****************************/ 
if  (shape  ==  1)  { 

Ipfconts  =  dvector(0,  numf irdelays) ; 

Ipfwts  =  dvector(0,  numf irdelays) ; 

ilpfconts  =  dvectord,  numfirtaps); 

ilpfwts  =  dvector(l,  numfirtaps); 

omega_c  =  setcutoff (samplerate,  lpfin_cutoff ) ; 

iomega_c  =  setcutoff (samplerate,  ilpfin_cutoff ) ; 

calcfilterweights  (numf irdelays,  omega_c,  Ipfconts,  Ipfwts); 
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calcf ilterweights  (numf irdelays ,  iojnega_c,  ilpfconts,  ilpfwts); 

} 

if  (carrierfreq  >  0.0)  { 

outweights  =  dvector(0,  numf irdelays) ; 

outcontents  =  dvectorCO,  numf irdelays) ; 

tiaf outweights  =  dvector(0,  numf irdelays) ; 

tiafoutcontents  =  dvector(0,  numf irdelays) ; 
nnwts  =  dvector(0,  numf irdelays) ; 
nnconts  =  dvector(0,  numf irdelays) ; 
omega_cl  -  setcutoff (samplerate,  lpfout_cutoff ) ; 
calcf ilterweights  (numf irdelays,  omega_cl,  outcontents, 
outweights) ; 

calcf ilterweights  (numf irdelays ,  oraega_cl,  tiafoutcontents, 
tiaf outweights) ; 

calcf ilterweights  (numf irdelays,  omega_cl,  nnconts,  nnwts); 

} 

tiafwts  =  dvector(0,  numtiafdelays) ; 

tiafconts  =  dvector(0,  numtiafdelays); 

mse  =  vector (0,  numsamples-l) ; 
tempmse  =  vector(0,  numsamples-l); 
tiafmse  =  vector(0,  numsamples-l); 
tiaftempmse  =  vector(0,  numsamples-l); 

for  (sample  =  0;  sample  <  numsamples;  ++sample)  { 
mse [sample]  =  0.0; 
tempmse [sample]  =  0.0; 
tiafmse [sample]  =  0.0; 
tiaftempmse [sample]  =  0.0; 

} 

inittiaf (tiafwts,  tiafconts,  numtiaftaps,  fetiaferror,  fetiafbias); 

tdafwts  =  inittdafwts(sampersym,  numtdaftaps) ; 

tdaf cents  =  inittdafeonts (numtdaftaps) ; 

tdaferror  =  initerror(sampersym) ; 

tdaf bias  =  initbias(sampersym) ; 

time  =  i£xnl(&nseed)  ; 

lasttime  =  time-datarate; 

itime  =  ran) (ftnseed) ; 

ilasttime  =  itime-idatarate; 

♦!»!****#*****  j|c**j(t**^!***!(t+*;(t********:|i*!)C)t:***!(e***+:****>(c  :*:**, *********** 

Simulation  Section 


c-i: 


4:!)c******:*!|c**:(t  !(!>(:  Jt:*************!!:**  **************  ***!)'=(!******  ****>!<*****!(:>(:**/ 


for  (loopnum  =  1;  loopnum  <=  nuinloops;  ++loopnum){ 
for  (sample  =  0;  sample  <  numsamples;  ++sample)  { 

/*********************************************************************** 
Interference  (SNOI)  Section 

***********************************************************************/ 
if  (icarrierampl  >  0.0)  { 

dlt  =  datagenCitime,  idatarate,  dlt,  feilasttirae, 

-  &bseed,  outputf lag) ; 

if  (format  ==  0) 

d2t  =  mai;'-hoster(dlt,  itime,  idatarate,  ilasttime); 
else 

d2t  =  bip  .»lar(dlt)  ; 
if  (shape  ==  1) 

sit  =  lpf(d2t,  numfirdelays,  ilpfuts,  ilpfconts,  igainl); 
else 

sit  =  d2t; 

if  (carrierfreq  >  0.0) 

s2t  =  modulate(slt,  icarrierfreq,  icarrierampl,  itime); 
else 

s2t  =  sit; 

} 

else 

s2t  =  0.0; 

/***********#*********************************************************** 

Signal  (SOI)  Section 

***************************+*+*****************************************/ 

d3t  =  datagen(time,  datarate,  d3t,  ftlasttirae, 
febseed,  outputf lag) ; 

if  (format  ==  0) 

d4t  =  mcinchester(d3t,  time,  datarate,  lasttime); 
else 

d4t  =  bipolar(d3t) ; 
s3tlast  =  s3t; 
if  (shape  ==  1) 

s3t  =  lpf(d4t,  numfirdelays,  Ipfwts,  Ipfconts,  gainl) ; 
else 
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1 


s3t  =  d4t; , 

I 

if  (carrierfreq  >0.0) 

s4t  =  modulate(s3t,  carrierfreq,  carrierampl ,  time); 
else 

s4t  =  s3t; 

Channel  Section 

+  **********************  ********=!!*♦***  *♦**!)=*****/ 
s5t  =  s4t  +  s2t  +  gasdev(&nseed)*noisegain; 

TDAF  Receix'er  Section 

)|c3t!***:|t:(c  +  !(c**********!t;**:j<*:t:  +  *i(t****!|c*-*******:|(j|t**%:jij(:****;(;:(c;(::jc!|cj(r+it:*  !*•**:(<:(<*:****  / 

s6t  =  tdafCsSt,  s4t,  mu,  sampersym,  numtdaftaps,  tdafwts, 
tdafconts,  tdaf error,  tdafbias,  frozen); 

I 

if  (carrierfreq  >  C) 

s7t  =  demodulateC  -t,  carrierfreq,  deltaphase,  time); 
s8t  =  lpf(s7t,  numfirdelays,  outweights,  outcontents, 
gain2) ; 

} 

else 

s8t  =  s6t: 

y  ic :te  :f(  %  :tc  3(c  9(c  ^  .ic  ;jc  :|e ^  3):  9|c  :fc  3k  ^  :4c  3^  sfc  ;4c  )fc  :ic :(( ;fc  s(c 

TIAF  Receiver  Section 

x6t  =  tiafCsSt,  s4t,  mu,  numtiaftaps,  &tiaf error, 
tiafwts,  tiafconts,  fetiafbias,  frozen); 

if  (carrierfreq  >  0.0)  { 

x7t  =  demodulate(x6t,  carrierfreq,  deltaphase,  time); 
x8t  =  lpf(}:7t,  numfirdelays,  tiaf outweights, 
tiaf outcontents,  gain2) ; 

} 

else 

x8t  =  x6t: 

I  :(c****!(c**!i!********4:***)!i****!(t!K****+**»****^J('*!!c********!tr!it*+***».***^!******^: 

Moise  Free  Section 

#********>!<****♦*******♦******♦********♦  >|t**!(:*****!|t**- *!«****!?  / 
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nn2tlast  =  nn2t; 
if  (carrierfreq  >  0.0)  ■[ 

r.nlt  =  demodulate(s4t,  carrierfreq,  deltaphase,  time); 
nn2t  =  lpf(nnlt,  numfirdelays,  nnHts,  nnconts,  gain2) ; 

} 

else 

nn2t  =  s4t; 

time  +=  samplerate; 
itime  +=  samplerate; 

/*:):**  Jit************  ***:(!!(!  +  *****♦:*;*********:)£**  *****  +  :([  :(c***!|c*;(c*:)t^e**** 

End  of  Receiver  section 

if  ((ioffsetfound)  &&  (loopnum  ==  adaploops  +!)){ 
if  (((nn2t  >  0.0)  &fe  (nn2tlast  <  0.0))  II 
((nn2t  <  0.0)  ft&  (nn2tlast  >  0.0)))  { 
bitvar  =  0; 
tdafbitnum  =  0; 
tiafbitnum  =0; 
dsrdbitnum  =  0; 
offsetfound  =  1; 

} 

} 

if  (loopnum  ==  nuraloops)  { 

fprintf  (dsrdvolt,  "'/,d  7,f\n",  sample,  nn2t) ; 

.  fprintf (tdafvolt ,  '"/.i  '/,f\n",  sa-mple,  s8t); 

fprintf (tiaf volt,  ‘"/.i  7.f\n",  sample,  x8t); 

} 

if  (loopnum  >  adaploops)  { 
frozen  =  1; 

tempmse  fsample] =sqr (s4t-s6t ) / desiredpower ; 
tiaf tempmse [sample] =sqr(s4t-x6t) /desiredpower ; 
if  (format)  { 

++bitvar; 

tdafbitnum  +=  s8t; 
tiafbitnum  +=  x8t; 
dsrdbitnum  +=  nn2t; 
if  (!  (bitvar  7,  sampersym))  { 

if  (((tdafbitnum  >  0.0)  &&  (dsrdbitnum  <  0.0))  II 
((tdafbitnum  <  0.0)  &&  (dsrdbitnum  >  0.0))) 
++tdafbiterror; 


if  (((tiafbitnum  >  0.0)  &&  (dsrdbitnum  <  0.0))  j I 
((tiafbitnum  <  0.0)  &&  (dsrdbitnum  >  0.0))) 
++tiafbiterror; 

++bitcount ; 
tdafbitnum  =  0.0; 
tiafbitnum  =0.0; 
dsrdbitnum  =  0.0; 

} 

}  /*  end  if (format)  */ 

}  /*  end  if (loopnum>adaploops)  */ 

}  /*  Ends  inner  FOR  loop  */ 

if  (loopnum  >  adaploops)  { 

++msecount ; 

for  (samplel  =  0;  sa.mplel  <  numsamples;  ++samplel)  { 
mseCsamplel]+=tempmse[samplel3 : 
tiafmse [samplel] +=tiaf 1 empmse [samplel] ; 

} 

}  /*  Ends  outter  FOR  loop  *! 

for  (sa...ple  =  0;  sample  <  numsamples;  ++sample)  { 

fprintf  (tdafvec,"’/,i  */.f\n",  sample,  (float)  mse[sample]/msecount) ; 
fprintf (tiafvec,"*/,i  '/.f\n",  sample, 

(float)  tiafmse [sample] /msecount  ); 

} 

avg  =  0.0; 
tiafavg  =  0.0; 
count  =  0; 

for  (sample  =  0;  sample  <  numsam.ples;  ++sample)  { 
avg  +=  (double)  mse [sample] ./msecount; 
tiafavg  +=  (double)  tiafmse [sample] /msecount; 

++count ; 

} 

avg/ =count ; 
t i af  avg/ =count; 

fprintf  (numbers  ."average  error  for  TDAF:  ’/,f\n",  avg); 
fprintf  (numbers,  "average  error  for  TIAF:  ’/.fXn",  tiafavg); 
fprintf  (numbers,  "number  of  tdaf  bit  errors:  */,d  ",  tdafbiterror); 
fprintf  (numbers,  "for  a  bit  error  rate  of  7.f\:i",  (float) 
tdafbiterror/bitcount) ; 

fprintf  (numbers,  "number  of  tiaf  bit  errors:  ’/,d  ",  tiaf  biterror) ; 
fprintf  (numbers,  "for  a  bit  error  rate  of  '/,f\n",  (float) 
tiaf biterror/bitcou It) ; 
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fclose(tdafvolt) ; 
fclose(t;afvolt) ; 
f close(dsrdvolt) ; 
f close (tdafvec) ; 
fclose(tiafvec) ; 
fclose(nuinbers) ; 

return  0; 

}  /*  Ends  mainO  */ 

/**s»;*******>t:*************!t:********5(!************************tf/ 

double  setcutoff (double  samplerate,  double  cutoff)  { 

double  temp  =  (2  *  pi  *  cutoff  *  samplerate) ; 
return  temp; 

} 

/***♦*>*!:(£:♦:!(:********!)£*  *#*)(!  +  **!(!★*****!('♦***»>)'*******♦****♦**♦***/ 

void  calcfilterweightsCint  numdelays,  double  omega_c, 
double  ♦contents,  double  ♦weights)  { 

int  count  =  -1; 

double  normf actor  =  1.0; 

double  temp,  window; 

double  ♦contents_ptr  =  contents; 

double  ♦weight„ptr  =  weights; 

double  *array_end  =  weights  +  nximdelays  +  1; 

double  H  =  numdelays; 

for  (  weight _ptr  =  weights;  weight_ptr  <  array_end;  ++weight_ptr  )  { 

count  +=  1; 

*eontents_ptr  =  1.0; 

++contents_ptr; 
temp  =  count  -  (H/2) ; 

window  =  0.54  -  0.46^cos(2*pi*count/(H)); 
if  (temp  ==  0)  { 

*weight_ptr  =  window^normf actor; 

} 

else 
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} 


*weight_ptr  =  window* (sin(oraega_c*temp))/(pi*temp)*normf actor; 


} 


#define  IBl  1 
#define  IB2  2 
#define  IBS  16 
ttdefine  IB18  131072 
#define  MASK  IB1+IB2+IB5 

int  irbit2(unsighed  long  *bseed)  { 

if  (*bseed  &  IB18)  { 

*bseed  =  ((  *bseed  "  MASK)  «  1)  I  IBl; 
return  1 ; 

} 

else  { 

*bseed  <<=  1; 
return  0; 

} 

} 


int  datagen (double  time,  double  datarate,  int  lastdata,  double 

*lasttime,  unsigned  long  *bseed,  int  outputflag)  { 


int  temp; 

if  (time  >=  (*lasttime  +  datarate))  { 
*lasttime  =  time; 
if  (outputflag  ==  1)  { 
temp  =  irbit2(bseed) ; 

} 

else 

temp  =  -lastdata; 

} 

else 
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} 


temp  =  lastdata; 
return  temp; 


double  Ipf (double  input,  int  numdelays,  double  ♦weights, 
double  *contents,  double  gain)  { 

int  count; 
double  sum  =  0; 

double  *contents_'^nd  =  contents  +  (numdelays); 
double  *weights_end  =  weights  +  (numdelays) ; 

for  (count  =  numdelays;  count  >  0;  — ^^count)  { 

♦(contents_end)  =  ♦(contents_end  -  1); 
sum  +=  *weights_end  ♦  * contents. end; 

— contents.end; 

— weights. end; 

} 

♦contents  =  input; 

sum  +=  ♦weights  ♦  ♦contents; 

return  sum  ♦  gain; 


} 


double  modulate (double  input,  double  carrierfreq, 
double  carrierampl,  double  time)  { 


double  temp  =  input+carrierampl^cos(2^pi^carrierfreq+time) 
return  temp; 


} 

/  *♦♦♦♦*♦**♦♦♦♦*;(■**♦♦♦♦♦**♦♦  =♦♦♦:>(!**/ 
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double  demodulate (double  input,  double  carrierfreq^ 
double  phase,  double  time)  { 

double  temp  =  2  *  input  *  cos (2  *  pi  *  carrierfreq  * 
time  +  phase) ; 

return  temp; 


} 

/*  The  following  routine  was  taken  from 
"Numerical  Recipes  in  C"  by  Press  et  al.  */ 

#define  Ml  259200 
#define  lAl  7141 
#define  ICl  54773 
#define  RMl  (1.0/Ml) 

#define  M2  134456 
#define  IA2  8121 
#define  IC2  28411 
#define  RM2  (1.0/M2) 
ftdefine  M3  243000 
#define  IA3  4561 
#define  IC3  51349 

float  ranKint  *idum)  { 


static  long  ixl,ix2,ix3; 
static  float  rC98]; 
float  temp; 
static  int  iff=0; 
int  j  ; 

void  nrerrorO; 

if  (*idum  <011  iff  ==  0)  { 
iff=l ; 

ixl=(ICl-(*idum))  '/,  Ml; 
ixl=(IAl>t=ixl+ICl)  y.  Ml; 
ix2=ixl  '/.  M2; 
ixl=(IAl*ixl+ICl)  y.  Ml; 
ix3=ixl  y,  H3; 
for  (j=l; j<=97;j++)  { 

ixl=(IAl=*=ixl+ICl)  y.  Ml; 
ix2=(IA2*ix2+IC2)  '/.  M2; 
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rCj]=(ixl+ix2*RM2)*RMl ; 

} 

*idum=l; 

} 

ixl=(IAl*ixl+ICl)  •/.  Ml; 
ix2=(IA2*ix2+IC2)  '/.  M2; 
ix3=(IA3*ix3+IC3)  */.  H3; 
j=l  +  ((97=!-.ix3)/M3); 

if  (j  >  97  I  1  j  <  1)  hrerrorC'RANl :  This  cannot  happen."); 
cemp=r[j] ; 

rCj]=(ixl+ix2*RM2)*RMl ; 
return  temp; 


#undef  HI 
Sundef  lAl 
#undef  ICl 
Sundef  RMl 
#undel  M2 
#undef  IA2 
Sundef  IC2 
#undef  RH2 
#undef  H3 
#undef  IA3 
#undef  ICS 


/*  The  following  routine  was  taken  from 
"Numerical  Recipes  in  C"  by  Press  et  al.  */ 

float  gasdev(int  +idum)  { 

static  int  iset=0; 
static  float  gset; 
float  fac,r,vl,v2; 
float  ranlO; 

if  (iset  ==  0)  { 
do  { 

vl=2.0*ranl(idum)-1.0; 
v2=2 . 0*ranl ( idum) -1.0; 
r=vl*vl+v2*v2; 

}  while  (r  >=  1.0  | |  r  ==  0.0); 
fac=sqrt(-2.0*log(r)/r) ; 
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gset=vl*fac; 

iset=l; 

return  v2*fac: 

}  else  { 
iset=0; 
return  gset; 

> 

} 

double  **inittdafuts(int  sampersym,  int  numtaps)  { 

int  rownum,  colnum; 
double  **tdafwts; 

tdafwts  =  dinatrix(l,  sampersym,  1,  numtaps); 
for  (rownum  =  1;  rownum  <=  sampersym;  ++rownum) 

for  (colnum  =  1;  colnum  <=  numtaps;  ++colnum) 
tdafwts [rownum]  [colnum]  =  0.0; 


return  tdafwts; 


} 

double  *inittdafconts(int  numtaps)  { 

int  tapnum; 
double  4:tdafconts; 
tdafconts  =  dvector(l,  numtaps); 
for  (tapnum  =  1;  tapnum  <=  numtaps;  ++tapnum) 
tdafconts [tapnum]  =  0.0; 
return  tdafconts; 

} 

dou.^le  tdaf (double  input,  double  desired,  double  mu, 
int  sampersym,  int  numtaps,  double  **tdafwts, 
double  *tdaf cents ,  double  *errbr,  double  *bias,  int  frozen)  { 
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input : 

desired: 

mu: 

sampersym 


/  This  routine  is  the  Time  Dependent  Adaptive  Filter. 

/ 

/  Input  Parameters'; 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 


The  value  to  be  filtered. 

The  training  signal. 

Used  to  calculate  the  filter  weights. 
The  number  of  samples  per  symbol  of 
the  input  baseband  data  sequence. 

The  TDAF  is  made  up  of  sampersym 
TIAFs . 

numtaps:  The  number  of  taps  in  each  of  the 

member  TIAFs  making  up  the  TDAF. 
tdaf wts :  An  n  by  m  array  where  n  =  numtaps  and 

m  =  sampersym  containing  the  current 
filter  coefficients  of  the  TDAF. 
tdaf cents:  A  vector  that  holds  the  past  values 
of  the  input. 

error:  A  vector  of  length 

sampersym  that  Contains  the  quantity 
"desired  -  sum"  (see  sum  below)  for  each 
TIAF  in  the  TDAF. 

Local  Variables: 

callnum:  Used  to  calculate  which  TIAF  supplies 

the  output  for  each  call  to  the  function 
tapnum:  Indexing  variable  for  the  taps  of  each 

TIAF. 

filternum:  The  index  that  points  to  the  TIAF 

that  supplies  the  output  of  the  function, 
sum:  Holds  the  intermediate  cind  final  value 

of  the  output  of  the  function. 


static  int  callnum  =  0; 
int  tapnum,  filternum; 
double  sum  =  0.0; 

filternum  =  (callnum  */,  sampersym)  +  1; 

++callnixin; 

/*  Shift  data  in  filter  to  the  right _ */ 

for  (tapnum  =  numtaps;  tapnum  >  1;  — tapnum)  { 
tdaf conts [tapnum]  =  tdaf conts [tapnum  -  1]; 

/*  and  calculate  the  output  due  to  the  past  values  */ 
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sum- +=  tdafcontsCtapnum]  *  tdafwts[filtern™3  [tapnum];; 

} 

/*  Now=  shift  the  input  into  tap  number  1  */ 
tdafcontsCl]  =  input; 

sum  +=  tdafcontsCl]  *  tdaf wts [f iltefnum]  [1] ; 

/*  Add  the  contribution  due  to  the  bias  weight  */ 
sum  +=  biasCfilternum] ; 

/*  Calculate  the  error  */ 
error [filternum]  =  desired  -  sum; 

/*  Update  the  f ilterweights  for  the  next  call  */ 
if  (Ifrozen)  { 

for  (tapnum  =  1;  tapnum  <=  numtaps;  ++tapnum)  { 

tdafwtsCf ilternum] [tapnum]  +=  2  *  mu  *  error [filternum] 
*  tdafcontsCtapnum]; 

} 

/*  Update  the  bias  weight  */ 
biasCfilternum]  +=  2  *  mu  *  error [filternum] ; 

} 

return  sum; 


y  ^  *  4: 4;  4:  4:  *  4:  If:  4:  ^  4= ’f:  *  4:  if:  4=  4:  *  4: 4:  *  4:  ^  ^  t  ^  ^  ^  ^  ^  ^  it:  / 

double  *initerror(int  sampersym)  { 

int  filternum; 
double  *error; 

error  =  dvectorCl,  sampersym); 

for  (filternum  =  1;  filternum  <=  sampersym;  ++f ilternum) 
error [filternum]  =  0.0; 
return  error; 

} 


/  4c4:4:4t4c4:4c4t4c4t4t4c4c4t4:4:44‘=(‘4‘4:4:4:4:4:4:4:4'4:4!4:4:4:4c4:*4'4'4:4'4:4t4:4:4:4:4'4:4t4:4:4'4:**4'4!4:4'/ 

double  *initbias(int  sampersym)  { 

int  filternum; 
dcuble  *bias; 

bias  ••=  dvectorCl,  sampersym); 

for  (f iiternum=l ;filternum<=sampersyiu;++f ilternum) 
biasCfilternum]  =  0.0; 
return  bias; 
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} 


/j(!*****!(c**!(:**:(c*!tt*;t!*^  •  ■i:*:*j|t***!|cj(:***********5)f******************/ 

/*  The  following  routine  was  taken  from 
"Numerical  Recipes  in  C"  by  Press  et  al.  */ 

void  nrerror(chaf  errbr_text [] )  { 

fprintf(stderr, "Numerical  Recipes  run-time  error. . An") ; 
f  printf  (stderr ,  "'/,s\n" ,  error  .text) ; 
fprintf (stderr," . . .now  exiting  to  system. . .\n") ; 
exit(l) ; 


/*  The  following  routine  was  taken  from 
"Numerical  Recipes  in  C"  by  Press  et  al.  */ 

double  *dvector(int  nl,  int  nh)  { 

double  *v; 

v=(double  *)malloc( (unsigned)  (nh-nl+l)*si2eof (double)) ; 
if  (!v.)  nrerrorC allocation  fa^ure  in  dvectorO"); 
return  v-nl ; 

} 

/j|c^tjf!*  +  *:):*j(!*******  +  ***^c**!|'******!t:*+******=l<*****!|=*******>l'  +  ****/ 

/*  The  following  routine  was  taken  from 
"Numerical  Recipes  in  C"  by  Press  et  al.  */ 

double  **dmatrix(int  nrl,  int  nrh,  int  ncl,  int  nch)  { 

int  i; 
double  **m; 

m=(double  **)  malloc( (unsigned)  (nrh-nrl+3.)+sizeof  (double*)) ; 
if  (!m)  nrerrorC'allocation  failure  1  in  dmatrixO"); 
m  -=  nrl; 

for(i=nrl;i<=nrh;i++)  { 

m[i]=(double  *)  malloc( (unsigned)  (nch-ncl+l)*sizeof (double)) ; 
if  (!m[i3)  nrerrorC'allocation  failure  2  in  dmatrixO"); 
m[i3  -=  ncl; 
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tempi  =  abs(input-l) ; 
else 


tempi  =  input; 
if  (tempi  ==  1) 
temp2  =  1; 
else 

temp2  =  -1; 
return  temp2; 

} 

/*  The  following  routine  was  taken  from 
"Numerical  Recipes  in  C"  by  Press  et  al.  */ 

void  free_dmatrix (double  **m,  int  nrl,  int  nrh,  int  ncl,  int  nch) 

int  i ; 

for(i=nrh;i>=nrl;i — )  free((char*)  (m[i]+ncl)); 
f ree( (char*)  (m+nrl) ) ; 

} 

/  j|c!t;i)!:)c3|t*l!*******!(tj(t****j|e*****j(c*3):5|c*)|:*j|t*!(e*****j|c***j)c***j)t!)tj|t!tc*)(t*!t;**  J 

/*  The  following  routine  was  taken  from 
"Numerical  Recipes  in  C"  by  Press  et  al.  */ 

void  free_dvector (double  *v,  int  nl,  int  nh) 

{ 

free((char*)  (v+nl)); 

} 

/:<c4:;<(**4!***************************************>!:*«!|:**********/ 

/*  The  following  routine  was  taken  from 
"Numerical  Recipes  in  C"  by  Press  et  al.  *f 
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i,ii„iiii  |i|i,;i,|i  |||'||ii'l|lill!lllll^|lJ|l!l||,|l|||||1llll:ll|||ili;llJ;il^llJ^  . . . .  I'liHiii  i:ii'Mi!ii"i:iii:iiiiiiiiiii!i  i:ii  lllllll‘■llill'll  . . .  !iii”UiiiL-ailiiiiii  iijia<aFh>  uii  iii 


float  *vector(int  nl,int  nh)  { 


float  *v; 

v=(float  *)malloc( (unsigned)  (nh-nl+l)*sizeof (floai 
if  (!v)  nrerror("allocation  failure  in  vectorO"); 
return  v-nl; 

} 

/****:tc*5t:******!(!**************:t:****i|:*3|t*!(:*+*J|c*)(c*!(t******  +  *!(c**j|c*/ 

void  inittiaf  (double  ^weights,  double  *contents,  int  numtaps, 
double  terror,  double  *bias)  { 


int  count; 

♦error  =0.0; 

♦bias  =  0.0; 

for  (count  =  0;  count  <  numtaps;  ++count)  { 

♦(weights  +  count)  =  0; 

♦(contents  +  count)  =  0; 

} 

} 

/*****  *!(<*^c****^**:jt****^:*  **:(!♦!(;******>(:  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦/ 

double  tiaf  (float  input,  float  desired,  float  mu,  int  numtaps, 

double  ♦error,  double  ♦weights,  double 
♦contents,  double  ♦bias,  int  frozen)  { 

int  numdelays  =  numtaps  -  1 ; 
int  count; 
double  sum  =0.0; 

for  (count  =  numdelays;  count  >  0;  — count)  { 

contents  [count]  =  contents [count  -  1]  ; 
sum  +=  contents  [count]  ♦  weights [count] ; 

} 


contents [0]  =  input; 

sum  +=  contents [0]  ♦  weights [0] ; 

sum  +=  ♦bias; 
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♦error  =  desired  -  sum; 


if  ('.frozen)  { 

for  (count  =  0;  count  <  numtaps;  ++count) 

weights [count]  +=  2  *  mu  *  *errof  ♦  contents [count] ; 
♦bias  +=  2  ♦  mu  ♦  ♦error; 

} 

xeturn  sum; 


} 

/  4:  ^  :tc  !|c  *  4:  :|c  :|c  :(c :)!  :jc  4: 4:  *  4: 4: 4:  ^  4:  :(<  ♦  4=  *  >1' *  4: 4: ’ll  4=  ^  *  :4c  *  *  4:  ♦  ♦  ♦  4:  *  / 

int  bipolar(int  input)  { 


if  (  input  ==  1  ) 
return  1; 
else 

return  -1; 

} 
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Appendix  D.  Source  Code  for  the  LC  Version 


/****  This  program  is  a  simulated  digital  communication  system  *******1 

/:tc^:|c:K*:4;**:t;4:;tc«:t!!ic^:t;:|c:|c4;i|ci|:*:t!*:t;:t:*****«**************«^************:ic4:!t:*«*ii:>l  */ 

/*  30  Sep  91:  Hade  number  of  symbols  for  simulation  user  selectable  / 
/  Executable  file:  Imcrv  / 

/**«*f:ic:ic*:)c:t;«4c:t;:tc:t:****«********:|‘4:***«***«**«***«*****«1c«!l‘*****^1‘ «-«******/ 

#define  spare  1 

#include  <stdio.h> 
ftinclude  <math.h> 

#include  <stdlib.h> 


^define  pi 
Sdefine  data.freq 
^define  data.rate 
ifdefine  sqr(x) 


(4*atanCl,0)) 

1.0 

1.0/datafreq 

(x)*(x) 


/■^*it*******************************************************/ 


/*  Function  Prototypes  */ 

/*  Refer  to  Appendix  A  -  All  Functions  are  the  same  */ 

/ifiiiHf**1f***ilfL:it:^tiif*^^fAf*******^******^^iii************************/ 

int  raainO  { 


unsigned  long  numsamples; 

int  numfirtaps,  numtiaftaps,  nuratdaftaps; 

int  numf irdelays ,  numtiaf delays ; 

int  numloops ; 

int  samplefreq; 

int  count ; 

int  outputflag; 

int  numsym; 

float  loopf actor; 

double  samplerate; 
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double  datafreq  =  data.freq; 

double  datarate  =  data_rate; 

double  carrierfreq; 

double  carrierampl ; 

double  noisegain; 

double  misadjust; 

double  mu; 

double  deltaphase; 

double  gainl; 

double  gain2; 

double  tiafgain2; 

duuole  avg  =  0.0; 

double  tiafavg  =  0.0; 

double  lpfin_cutoff ,  lpfout_cutoff ,  ilpf in_cutoff ; 
char  buffer [128]; 

double  *lpfconts,  *lpfHts; 

double  *nnwts,  *nnconts; 

double  *tiafwts,  *tiafconts; 

double  *outweights,  *outcontents; 

double  *tiafoutweights,  *tiafoutcontents; 

double  *mse,  *tempmse; 

double  *tiafmse,  *tiaftempmse; 

double  tiaf error  =  0.0; 

double  tiaf bias  =  0.0; 

int  sample,  samplel,  loopnum; 

int  dlt=l,  d2t,  dSt*!,  d4t,  dSt,  d6t; 

int  sampersym; 

int  format; 

int  shape; 

double  omega_c; 

double  omega_cl; 

double  time  =  0.0; 

double  sit,  s2t,  s3t,  s4t,  s5t,  s6t,  s7t,  s8t; 
double  x6t,  x7t,  x8t; 
double  nnlt,  nn2t; 
unsigned  long  int  bseed; 

int  nseed; 

double  lasttime  =  -datarate; 

double  **tdaf-wts ; 
double  *tdafconts; 
double  *tdaf error; 
double  *tdafbias; 


double  *tdafsyinbolmse: 
double  *tiafsymbolmse; 
int  msecount; 
int  msesize; 
int  mseindex  =  0; 
int  numsymbols; 
double  idatarate; 
double  idatafreq; 
double  ilasttime; 
unsigned  long  int  ibseed; 
double  *ilpfconts,  *ilpfwts; 
double  igainl; 
double  icarrierfreq; 
double  icarrierampl; 
double  itime  =  0.25; 
double  iomega.c; 
float  adapf actor  =  0.1178; 
double  dosiredpower; 
double  interf power; 
int  loopcount  =  0; 

FILE  *tdaflrn,  *tiaflrn,  *lnumbers: 

if  ((tdaflvn  =  fopenC'tdaflrn.dat" ,  "w"))  ==  NULL) 
printf  ("  ***  Could  not  open  tdaflrn.dat!  ***  \n"): 

if  ((tiaflrn  =  fopen("tiaflrn.dat" ,  "w"))  ==  NULL) 
printf  ("  ***  Could  not  open  tiaflrn.dat!  *=♦=*  \n"); 

if  ((Inumbers  =  fopenC'lnumbers.tex",  "w"))  ==  NULL) 
printf  ("  Could  not  open  lnunibers.tex!  ***  \n"): 


printf ("Seed  the  random  bit  generator:  "); 
gets (buffer) ; 

sscanf  (buffer,  "'/,U",  &bseed) ; 
printf  ("'/,d\n" ,  bseed) ; 

printf ("Seed  the  AWGK  Generator  (integer  <  0):  "); 
gets (buffer) ; 

sscanf  (buffer,  "'/,d",  &nseed) ; 
printf  ("'/,d\n" ,  hseed) ; 

printf ("Number  of  samples  per  symbol  required:  "); 
gets (buffer) ; 
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sscanf  (buffer,  "’/.d"  ,  fesampersyra) ; 
printf  ("'/,d\n" ,  sampersym) ; 

samplefreq  =  sampersym*datafreq; 
samplerate  =  (double)  1 .0/samplefreq; 

printf ("Manchester  (0)  or  Bipolar  (1)  format:  "); 
gets (buffer) ; 

sscanf  (buffer,  '"/.d",  feformat); 
printf  ("'/,d\n" ,  format) ; 

printf ("Pulse  shaping?  (l=y,  0=n) :  "); 
gets (buffer) ; 

sscanf  (buffer ,  "*/,d",  &shape) ; 
printf  ("y,d\n" ,  shape) ; 

printf ("Number  of  taps  in  the  FIR  filters:  "); 
gets (buffer) ; 

sscanf  (buffer,  '"/.d",  ftnumf  irtaps)  ; 
printf  ("y,d\n" ,  numf  irtaps) ; 

numfirdelays  =  numf irtaps  -  1; 

printf ("Number  of  taps  in  the  TIAF  adaptive  filter:  "); 
gets (buffer) ; 

sscanf  (buffer,  '"/.d",  fenumtiaftaps) ; 
printf  ("y,d\n" ,  numtiaftaps) ; 

numtiafdelays  =  numtiaftaps-1 ; 

printf ("Number  of  taps  in  the  TDAF  adaptive  filter:  "); 
gets (buffer) ; 

sscanf  (buffer,  "'/.d",  ftnurntdaf  taps) ; 
printf  ("y,d\n" ,  numtdaftaps) ; 

printf ("SOI  pulse  shaping  LPF  cutoff  (Hz):  "); 
gets (buffer) ; 

sscanf  (buffer,  '"/.If",  &lpf  in.cutoff ) ; 
printf  ("y,f\n" ,  lpfin_cutoff ) ; 

ilpf in_cutoff  =  lpfin_cutoff ; 

printf ("Gain  of  pulse  shaping  LPF:  "); 
gets (buffer) ; 

sscanf  (buffer,  '"/.If",  fegaihl); 
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printf  ("7.f\n" ,  gainl) ; 
igainl  =  gainl; 

printf ("SOI  carrier  amplitude:  "); 
gets (buffer) ; 

sscanf  (buff  er,  "'/,lf",  &carrierampl) ; 
printf  ("7, f\n"  ,  carrierampl) ; 

desiredpower  =  sqr(carrierampl)/2; 

printf ("SOI  carrier  frequency:  "); 
gets (buffer) ; 

sscanf  (buff  er,  "7, If",  ftcarrierf  req) ; 
printf  ("7.f\n"  ,  carrierfreq) ; 

printf ("Symbol  frequency  for  the  intorferer 
gets (buffer) ; 

sscanf  (buffer,  "7.1f",  feidatafreq); 
printf  ("7.f\n"  ,idatafreq) ; 

idatarate  =  1 .0/idatafreq; 

printf ("SNOI  carrier  amplitude:  "); 
gets (buffer) ; 

sscanf  (buffer,  "7, If",  &icarrierampl) ; 
printf  ("7.fAn"  ,  icarrierampl) ; 

interfpower  =  sqr(icarrierampl)/2; 

printf ("SNOI  carrier  frequency:  "); 
gets (buffer) ; 

sscanf  (buff  er,  "7ilf",  feicarrierfreq) ; 
printf ("7, f\n",  icarrierfreq).; 

printf ("Gain  of  output  LPF:  ") ; 
gets(buffer); 

sscanf  (buff  er,  "7,lf",  again2);- 
printf  (“7.f  \n" ,  gain2) ; 

printf ("Output  LPF  cutoff  (Hz):  "); 
gets (buffer) : 

sscanf  (buffer.,  "7, If",  &lpf  out.cutoff  )  ; 
printf  ("7.f\h" ,  Ipfout^cutoff ) ; 


prihtf ("Phase  shift  for  demodulator:  "); 
gets (buffer) ; 

sscanf  (buffer,  "’/.If",  &deltaphase) ; 
prihtf  ("7, f\n" ,  deltaphase) ; 

prihtf ("Noise  factor:  "); 
gets (buffer) ; 

sscanf  (buffer ,  "7, If",  ftnoisegain) ; 
printf  ("7.f\n" ,  noisegain) ; 

prihtf ("Misadjustment  factor:  "); 
gets(buff er) ; 

sscanf  (buffer,  "7.1f",  femisadjust); 
printf  ("7.f\n" ,  misadjust) ; 

if  (carrierfreq  >  0.0) 

mu  =  misadjust/((desiredpower  +  interfpower  + 
sqr (noisegain) )*(numtdaf taps)) ; 

else 

mu  =  misadjust/((sqr(carrierampl)  +  sqr(icarrierampl)  + 
sqr (noisegain) )*(numtdaf taps)) ; 

printfC'Mu  =  7..10g\h",,  mu); 
f  printf  (Inumbers,  "Mu  =  7..10g\n",  mu); 

prihtf ("Random  data  (1)  or  square  wave  (0):  "); 
gets (buffer) ; 

sscanf  (buffer ,  "7.d"  ,  feoutputf  lag) ; 
printf  ("7.d\n" ,  outputf  lag) ; 

printf ("Number  of  epochs:  "); 
gets (buffer) ; 

sscanf  (buffer,  "7.d",  fenumloops); 
printf  ("7.d\n" ,  numloops); 

printf ("Number  of  symbols  to  average:  "); 
gets (buffer) ; 

sscanf  (buffer,  "7.d",  ftnumsym) ; 
printf  ("7.d\n"  ,  nurasym) ; 

printf ("Number  of  symbols  in  simulation:  "); 
gets (buffer) ; 

sscanf  (buffer,  "7.d",  ftnumsymbols)  ; 
prihtf  ("7id\n" ,  numsymbols) ; 
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numsamples  =  numsymbols  *  sampersym; 

fprintf (Inumbers, "Number  of  symbols  in  learing  curve:  */,d\n", 

numsymbols) ; 


loopfactor  =  1 .0/(numloops) ; 
if  (shape  ==  1)  { 

Ipfconts  =  dvector(0,  numf irdelays) ; 

Ipfwts  =  dvector(0,  numf irdelays) ; 

ilpfconts  =  dvectord,  numfirtaps); 

ilpfwts  =  dvectord,  numfirtaps): 

omega_c  =  setcutoff (samplerate,  Ipf in_cutoff ) ; 

iomega_c  =  setcutoff (samplerate,  ilpf in_cutoff ) ; 

calcf ilterweights  (numf irdelays,  omega_c,  Ipfconts,  Ipfwts); 

calcf iltefweights  (ninr.f irdelays,  iomega_c,  ilpfconts,  ilpfwts); 

} 

if  (carrierfreq  >  0.0)  { 

outweights  =  dvector(0,  numf irdelays) ; 

outcontents  =  dvector(0,  numf irdelays) ; 

tiafoutweights  =  dvector(0,  numf irdelays) ; 

tiaf outcontents  =  dvector(0,  numf irdelays) ; 

nnwts  =  dvector(0,  numf irdelays)  ; 
nnconts  =  dvector(0,  numf irdelays) ; 
omega.cl  =  setcutoff (samplerate,  Ipfout.cutoff ) ; 
calcf ilterweights  (nximf irdelays,  omega.cl,  outcontents, 
outweights) ; 

calcf ilterweights  (numf irdelays,  omega.cl,  tiaf outcontents , 
tiafoutweights) ; 

calcf ilterweights  (numf irdelays,  omega.cl,  nnconts,  nnwts); 

} 

tiafwts  =  dvector(0,  numtiafdelays) ; 

tiafconts  =  dvector(0,  numtiafdelays); 

msesize  =  numsamples/sampersym; 
mse  =  dvector(0,  msesize  -  1); 
tempmse  =  dvector(0,  msesize  -  1); 
tiafmse  =  dvector(0,  msesize  -  1); 
tiaf tempmse  =  dvector(0,  msesize  -  1); 
tdafsymbolmse  =  dvector(0,  sampersym-1) ; 
tiafsymbolmse  =  dvector(0,  sampersym-1); 

for  (sample  =  0;  sample  <  msesize;  ++sample)  { 
rase [sample]  =  6.0; 
tempmse [sample]  =  0.0; 


tiafmse [sample]  =  0.0; 
tiaftempmse[sOTiple]  =  0.0; 

} 

for  (loopnum  =  1;  loopnum  <=  numloops;  ++ioopnum){ 

inittiaf (tiafwts,  tiafconts,  numtiaftaps,  fttiaferror,  &tiafbias) ; 
tdafwts  =  inittdafwts(sampersym,  nvuntdaftaps) ; 
t.dafconts  =  ihittdafconts(numtdaftaps) ; 
tdaferror  =  initerror(sampersym) ; 

'dafbias  =  initbias(sainpersym) ; 
msecount  =  0; 
time  =  ranlC&nseed) ; 
lasttime  =  time-datarate; 
itime  =  ranl(&hseed) ; 
ilasttime  =  itime-idatarate ; 

for  (sample  =  0;  saiiiple  <  numsamples;  ++sample)  { 

Interference  (bNOI)  Section 

if  (icarrierampl  >  0.0)  { 

dlt  =  datagenC itime,  idatarate,  dlt,  feilasttime, 
ftibseed,  outputflag).; 


if  (format  =~  0) 

d2t  -  raanchester(dlt,  itime,  idatarate,  ilasttime); 
else 

d2t  =  bipolar(dlt) ; 
if  (shape  ==  1) 

sit  =  lpf(d2t,  numf irdelays ,  ilpfwts,  ilpfconts,  igainl); 
else 

sit  =  d2t; 

if  (carrierfreq  >  0.0) 

s2t  =  modulate(slt,  icarrierfreq,  icarrierampl,  itiasa) ; 
else 

s2t  =  sit; 

} 

else 

s2t  =  0.0; 

Signal  (SOI)  Section 
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d3t  =  detagen(time,  datarate,  d3t,  felasttime, 

&bseed,  outputf lag) ; 

if  (format  ==  0) 

d4t  =  manchester(d3t,  time,  datarate,  iasttimeV ■ 
else 

d4t  =  bipolar(d3t) ; 
if  (shape  ==  1) 

s3t  =  lpf(d4t,  iiumf irdelays,  Ipfwts,  Ipfconts,  gai;:i.> . 
else 

s3t  =  d4t; 

if  (carrierfreq  >  0.0) 

s4t  =  modulate(s3t ,  carrie-rf req ,  carrierampl ,  time); 
else 

s4t  “  £3t; 

^  ^  sf:  If.' ^  ^  ^  :fc  ^  :fc  3(c  :fc  ^  ^  :4e  4:  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 

Channel  Section 

3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3fC  3^  3^  3^  3fc  3^  3^  3^  3^  3^  3f(  3^^  3^  3^  3^  3^  3^  3fC  3|C  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  3^  J 

s5t  =  s4t  +  s2t  +  gi;sdev(&nseed)*noisegaj.n ; 

TDAF  Receiver  Section 

s6t  =  tdaf(.s5t,  s4t,  mu,  sampersym,  numtdaftaps,  tdafwts, 
tdafconts,  tdaf error,  tdafbias) ; 

if  (carrierfreq  >  0)  { 

s7t  =  demodulate (s6t,  carrierfreq,  deltaphase,  time); 
s8t  =  lpf(s7t,  numf irdelays,  outweights,  outcontents, 
gain2) ; 

} 

else 

s8t  =  s6t; 

/.  3t(  3f3  34(  34(  3f(  3f(  3f(  3|(  3);  3{C  3f(  3^  3|C  3f:  ^  3f(  3^  3^  3f:  3tc  3i(  3}(  3K  3f:  3t:  3f(  3fC  3^  3{C  S(C  3fC  3(C  3^  3i(  3fC  3|C  3((  3^  3)C  3(C  3|C  3fC  3t .  3f(  3t:  3f(  3fC  3|C  3**  ie  3(C  3*;  3);  3^  3^  3fe  .fC 

TIAF  Receiver  Section 

3f(3fc3|c:f(3f(34C3((3f(3f(34:3(:3|e3f;3t;3f(.'^3f(3fC3{(3^:|::f;3fe3K3te3|c3fC3^3t(3^3fe3(C3|(3i(3(e3{C3fC3|C3f(3tc3f:3^3(C3|C3fc:^3f(:fC3f(3fC3(C3tC3|(3t(3f:3f:)j;3f(i( ! 
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x6t  =  tiaf(s5t,  s4t,_mu,  numtiaf taps,  fttiaferror, 
tiafwts,  tiaf-bnts,  &tiaf bias) ; 

if  (carrierfreq  >0.6)  { 

x7t  -  demodulateCxOt,  carrierfreq,  deltaphass,  time); 
x8t  =  lpf(x7t,  numfirdeiays,  tiafoutweights, 
tiafoutcontents,  gain2) ; 

} 

else 

x8t  =  x6t; 

Noise  Free  Section 
if  (carrierfreq  >  0.0)  { 

nnlt  =  demodulate (s4t ,  carrierfreq,  deltaphase,  time); 
nn2t  =  lpf(nnlt,  numfirdeiays,  nnwts,  nnconts,  gain2) ; 

} 

else 

nn2t  =  s4t; 
time  +=  sampler ate; 
itime  +=  samplerate; 

mseindex  =  sample  7,  sampersym; 

tdafsymbolmseCmseindex]  =  (sqr(s4t-s6t))/(desiredpower) ; 
if  (mseindex  ==  (sampersym  -  1))  { 
avg  =0.0; 

for  (samplel  =  0;  sample!  <  sampersym;  ++samplel)  { 
avg  +=  tdafsymbolmse [samplel] ; 

} 

tempmse [rasecount]  =  avg/sampersym; 

++msecount; 

} 

if  (sample  <  msesize) 

tiaftempmse [sample]  =  (sqr(s4t  -  x6t))/(desii-::dpower) ; 

}  !*  Ends  inner  FOR  loop  *! 

++loopcount ; 

for  (samplel  =  0;  samplel  <  msesize;  ++samplel)  { 
mse [samplel] +=temprase [samplel] ; 
tiafmse [sample!] +=tiaf tempmse [samplel] ; 

} 

free^dmatrix(tdafwts ,  1,  sampersym,  1,  numtdaftaps) ; 
free^dvector(tdafconts,  1,  numtdaftaps); 
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free^dvector (tdaf error ,  1,  sampersym) ; 
frGe_dvector(tdafbias,  1,  sampersym); 

}  /*  Ends  outter  FOR  loop  */ 

for  (sample  =  0;  sample  <  msesize;  ++s^ple)  { 

fprintf(tdaflrn,"*/,i  */,f\n",  sample,  (float)  mse  [s  ample]// loop  count ) ; 
fprintf (tiaflfn,"’/,!  */.f\n",  sample, 

(float)  tiafmse [sample] *loopf actor) ; 

} 

avg  =  0.0; 
tiafavg  =  0.0; 
count  =  0; 
for 

(sample  =  0;  sample  <  msesize-numsym;  ^+sample)  { 
avg  +=  mse  [sample] /loopcount; 
tiafavg  +=  tiafmse  [sample] /loopcount; 

++count ; 

} 

avg/=count; 

tiafavg/=count; 

fprintf  (Inumbers,  "average  error  for  TDAF  over  ’/,d  samples:  '/,f\n", 
numsym,  avg); 

fprintf  (Inumbers,  "average  error  for  TIAF  over  */,d  samples:  '/,f\n", 
numsym,  tiafavg); 


fclose(tdaflrn) ; 
fclose(tiaflrn) ; 
fclosednumbars) ; 

return  0; 


}  /*  Ends  mainO  */ 
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